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ABSTRACT 

A pQCD-based model for parton production and equilibration in 
ultrarelativistic heavy-ion collisions is reviewed. The model combines 
pQCD processes including initial and final state radiations together 
with string phenomenology for nonperturbative soft processes. Nu- 
clear effects on the initial parton production, such as multiple parton 
scattering and nuclear shadowing of parton distribution functions are 
considered. Comparisons with existing data are made and further tests 
of the model to constrain model parameters are proposed. With the 
obtained space-time history of the parton production, evolution of the 
minijet gas toward a fully equilibrated parton plasma is studied. Direct 
probes of the early parton dynamics, such as pre-equilibrium photon 
and dilepton production, open charm production, J/ip suppression and 
jet quenching are also reviewed. 
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1 Introduction 



According to the Big Bang theory, during the early stage of the evolution of our 
universe, there existed a state of matter consisting mainly of unbound quarks and 
gluons at a temperature of about 200 MeV. As the universe expanded and cooled 
down, it went through a phase transition and eventually evolved to what we see today 
Such a phase transition was predicted fl], [2| shortly after Quantum Chromodynamics 
(QCD) was established as the fundamental theory for strong interactions. In this 
theory, quarks, which have three types of color charges, interact via non-Abelian 
gauge bosons, i.e., gluons. Because of the non-Abelian nature of strong interactions, 
quarks and gluons are confined to hadrons under normal circumstances. However, 
hadronic matter under extremely dense and hot conditions will go through a phase 
transition to form a quark-gluon plasma in which quarks and gluons are no longer 
confined to the size of a hadron. Lattice gauge studies of QCD at finite temperatures 
have indeed found such a phase transition, though the nature of the transition has 
not yet been clarified [[|. 

To recreate the early universe and experimentally verify the QCD phase tran- 
sition in the laboratory, it was proposed that nuclei be accelerated to extremely high 
energies and then allowed to collide with each other. During the collisions, both 
high baryon densities and high temperatures might be reached. Thus, a quark-gluon 
plasma could be formed which will subsequently go through the phase transition and 
hadronize into particles in the final states. During the last decade, heavy-ion experi- 
ments have been done at the AGS of BNL and the SPS of CERN with m 14 and 
200 AGeV, respectively. The results from those experiments ||] have demonstrated 
extremely rich physics which cannot be explained by simple extrapolation of pp col- 
lisions. However, up to now, there is no unambiguous evidence for the existence of a 
quark-gluon plasma over a significantly large space-time region. While experiments 
at AGS and SPS and their analyses continue, new experiments have been planned 
at RHIC of BNL and LHC of CERN with E cm = 200 AGeV and 5.5 ATeV, respec- 
tively. What I would like to emphasize in this review is that heavy-ion collisions at 
these ultrarelativistic energies will demonstrate a completely new dynamics which is 
not accessible at the present energies. 

The key issue here is the nuclear structure at different scales. When the 
transverse momentum transfer involved in each nucleon-nucleon collision is small 
(px < Aq CD ), effective models based on, e.g., meson-exchange and resonance forma- 
tion are sufficient to describe multiple interactions between hadrons in which the 
parton structure of the hadrons cannot yet be resolved. These coherent (with respect 
to the partons inside the hadron) interactions lead to collective behavior in low-energy 
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heavy-ion collisions as were first observed in the Bevalac experiments and recently 
at the AGS energies J/J. However, when px becomes large enough to resolve individ- 
ual partons inside a nucleon, the dynamics is best described on the parton level via 
perturbative QCD (pQCD). Though hard parton interactions occur at CERN-SPS 
energies (yfs ^ 20 AGeV), they play a negligible role in the global features of heavy- 
ion collisions. However, at collider energies (y/s ^ 100 AGeV) the importance of hard 
or semihard parton scatterings is clearly seen in high-energy pp and pp collisions f| . 
They are therefore also expected to be dominant in heavy-ion collisions at RHIC and 
LHC energies || These hard or semihard interactions happen on a very short 
time scale and they generally break color coherence inside the individual nucleons. 
After the fast partons pass through each other and leave the central region, a dense 
partonic system will be left behind which is not immediately in thermal and chemical 
equilibrium. The partons inside such a system will then further interact with each 
other and equilibration will eventually be established if the interactions are frequent 
enough among a sufficiently large number of initially produced partons. Due to the 
asymptotic behavior of QCD, production rates of hard and semihard partons are cal- 
culable via perturbative QCD during the initial stage of heavy-ion collisions. The 
color screening mechanism in the initially produced dense partonic system makes it 
also possible to use pQCD to investigate thermal and chemical equilibration of the 
system. 

Currently, there are many models to describe the coherent hadronic interaction 
in which parton degrees of freedom are not important yet. In these models, hadrons 
remain as individual entities and particle production is mainly through resonance 
formation. Indeed, most of the experimental data at AGS energies can be explained 
by RQMD |I| and ARC pf Monte Carlo simulations. As the colliding energy in- 
creases, partons inside nucleons become more relevant as the basic constituents of the 
interaction. However, nonperturbative soft interactions still dominate the collision 
dynamics in the energy range \fs ^ 100 AGeV. In order to describe particle produc- 
tion in this energy range, phenomenological string models |13|, 14, [15|, HJ have been 



developed and they can explain well the global properties of particle production at 
CERN SPS energies. The physics behind these string models is that multiple partons 
produced in soft interactions are not independent of each other. Coherence among 
them can be modeled by strings or flux tubes between leading quarks and diquarks (or 
antiquarks). Particle production comes from the fragmentation of the string through 
qq creation inside the flux tube. At collider energies, y/s > 100 AGeV, hard and 
semihard scatterings, which can be calculated via pQCD, become increasingly im- 
portant. Since nonperturbative physics is always present in any processes involving 
strong interactions, pQCD has to be combined with models of nonperturbative inter- 
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actions, e.g., string models. Such pQCD-inspired models, like PYTHIA JET], ISAJET 



IS], describe well the high phenomena in high-energy hadronic collisions. Along 



the same line, M. Gyulassy and I extended this approach to high-energy heavy-ion 
collisions and developed the HIJING (Heavy Ion Jet INtercation Generator) Monte 
Carlo model |L9| to take into account the physics of hard and semihard parton pro- 



duction. Nearly at the same time, similar models, PCM f2(| and DTUNUC [21], also 
emerged. In particular, PCM has explicitly modeled the space-time evolution of the 
parton production processes and rescatterings along the lines of Boal f22fl . Though 
hard and semihard processes become dominant at high energies, they are always ac- 
companied by soft interactions. These soft interactions in general produce many soft 
final partons and take longer time (~ 1 fm/c) to complete. Therefore, they cannot be 
modeled by simple elastic parton-parton scatterings. However, these nonperturbative 
interactions might be modified inside a dense partonic medium, so that a perturbative 
method might be sufficient to describe the evolution of the initially produced dense 
partonic system towards an equilibrated quark-gluon plasma. 

The crux of parton dynamics in heavy-ion collisions rests in the complexity of 
multiple parton scatterings and the associated bremsstrahlung. An exact quantum 
field treatment of these multiple parton interactions must include all possible matrix 
elements and the interference between them. Although quantum transport theory 

|23], |25[ can in principle include all these interference effects and developments in 
this subject are impressive in the last few years, its applicability to heavy-ion collisions 
is still far from being realistic. An alternative and more realistic solution to this 
problem is to find ways to incorporate correctly the subtle interference phenomena in 
a classical parton cascade model. There are two important interference phenomena 
in multiple parton interactions. (1) The interference among different amplitudes 
of multiple parton scatterings, especially in the initial parton scatterings when the 
two beams of fast partons pass through each other, leads to the Glauber formula 
of multiple interactions |26|, pTfl . Thus, a semiclassical cascade cannot be used to 
treat the initial parton production during the overlapping period of the two colliding 
nuclei. The resultant interference is responsible, as I will demonstrate later, for the 
disappearance of the nuclear enhancement of jet production in pA collisions at large 
Pt and high energies. Therefore, we can neglect processes in which a parton suffers 
multiple hard scatterings in the initial parton scatterings. Brodsky and Lu have shown 

that a Glauber analysis of multiple antiquark-nucleon scatterings also leads to 
nonperturbative shadowing of the quark distributions inside heavy nuclei. (2) The 
spectrum of induced gluon bremsstrahlung in multiple scatterings is also modified 
by the destructive interference among different radiation amplitudes, the so-called 
Landau- Pomeranchuk-Migdal (LPM) effect [EO, [50fl. In a detailed analysis W\ 
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one finds that the interference actually happens between two amplitudes in which 
the beam parton has completely different virtualities, i.e., time-like in the final state 
radiation of one scattering and space-like in the initial state radiation of the previous 
one. Therefore, a propagating parton can no longer be considered as always time-like 
with a decreasing virtuality in a parton cascade simulation, e.g., PCM J20| . 



Another aspect of parton production and interaction in ultrarelativistic heavy- 
ion collisions is that the same processes responsible for parton equilibration can also 
provide direct probes of the early parton dynamics and the properties of the parton 
gas during its evolution toward a fully equilibrated quark-gluon plasma. For exam- 
ple, parton scatterings in an equilibrating system also lead to dilepton, photon and 
charmed quark production. Since the production rates of these three processes are 
proportional to the product of quark-quark, quark-gluon and gluon-gluon densities, 
respectively, measurements of the thermal enhancement of dileptons, photons and 
charmed quarks can provide us with information of the quark and gluon densities 
in the plasma, and thus on the thermal and chemical equilibration. Interactions of 
an initially produced high-p^ jet with partons inside the equilibrating plasma, in the 
meantime, will cause the jet losing energy. Studying of jet quenching due to energy 
loss can thus tell us whether and how the the parton gas is thermalized. Because of 
the same mechanism, J/ip is also suppressed due to its interaction with the parton 
gas. Since the J/ip dissociation cross section inside a deconfined partons is very dif- 
ferent from that inside a hadronic gas, the study of J/i/j suppression can then provide 
us with evidence for color deconfinement in the parton plasma and possibly with 
information on the QCD phase transition. 

This review will be structured as follows: In the next section, I will review a 
pQCD-inspired model for multiple parton production in pp, pA and AA collisions. 
Special emphasis will be devoted to the role of hard and semihard parton collisions 
and their connection to Pomeron structure and soft interactions. I will also discuss 
the interpolation between hard and soft physics and the possible treatment of non- 
perturbative interactions. In Section 3, nuclear effects such as multiple initial and 
final state interactions and nuclear shadowing of parton distribution functions will 
be discussed. The Monte Carlo implementation of parton production will also be 
reviewed. In Section 4, I will discuss the evolution of a partonic system via pQCD, 
including the Landau-Pomeranchuk-Migdal (LPM) effect on induced gluon radiation. 
In Section 5, I will discuss the possibility of using hard processes as probes of the 
early parton dynamics, the formation of quark-gluon plasma formation and the QCD 
phase transition. In particular, photon and dilepton production, charm production, 
pre-equilibrium J/ijj suppression, jet quenching, and monojet production will be dis- 
cussed. Finally I will give a summary and discussion in Section 6. 
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2 Parton Scatterings in Hadron-Hadron Collisions 



High-energy nuclear collisions are expected to be dominated by semihard parton col- 
lisions with transverse momentum transfers Pt ~ Po ~ 1 — 2 GeV/c. The produced 
partons carrying such transverse momenta are often referred to as minijets. Multiple 



parton production has been estimated [10| to produce up to 50% (80%) of the trans- 



verse energy per unit rapidity in the collisions of heavy nuclei at RHIC (LHC) energies. 
While not resolvable as distinct jets, minijets are expected to lead to a wide variety of 
correlations among observables, such as transverse momentum, strangeness and fluc- 
tuation enhancements, that compete with expected signatures of a QGP. Therefore, 
it is especially important to calculate these background processes as reliably as pos- 
sible. In addition, it has been shown that multiple mini-jet production is important 
in pp interactions to account for the increase of the total cross section^] with en- 
ergy, the increase of the average transverse momentum with charged multiplicity ||34|l . 
and the violation of Koba-Nielsen-Olesen (KNO) scaling of the charged multiplicity 



distributions |T7|, |35 



It has long been recognized that the inclusive jet cross section <7j e t in pp 
or pp collisions increases very rapidly with energy and eventually will be larger than 
the total inelastic cross section <Ji n . This is because the inclusive jet cross section also 
contains the average number of jets per inelastic event. Thus, 

CTjet = (njet)Cin, (1) 

where (^j e t) is the average number of hard or semihard scatterings per inelastic event. 
One way to include hard or semihard processes in the unitarized cross sections of pp 
and pp is to introduce the idea of Pomeron exchanges in an eikonal formalism. I will 
argue that hard or semihard scatterings can be considered as the hard loops inside 
a Pomeron exchange, and the interpolation between hard and soft components will 
become clear. This will also shed some light on the modeling of the soft component. 

2.1 Pomeron Exchange and Minijet Production 

The Pomeron was introduced to describe the effective interaction between quarks 



and gluons |36|, [37|, HH- ft nas been used to model hadronic cross sections, diffractive 
interactions and multiple particle production. In general, it is considered to be a 
color singlet object with the quantum numbers of a photon. It can couple both to 
quarks and gluons. The imaginary part of a Pomeron exchange amplitude is usually 
related to multiple particle production in hadronic interactions. Let us first consider 
the eikonal formalism for quark scatterings via multiple Pomeron exchange. 
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Denote the matrix element for a single Pomeron exchange as 

= 2m5(E l - E f )2nyfif x {s,t), (2) 



= ^^T( S ,q), (3) 
2ir iyj s 

where y/s/2 = E{ = Ef and g p is a coupling constant between a quark and a Pomeron. 
The amplitude f(s,t) is defined such that the differential cross section for elastic 
scatterings is given by 

^ = vr|/(,,t)| 2 . (4) 

The initial and final polarizations, <jj, <jf, should be averaged and summed over in 
the calculation of the cross sections. One can check that if A(q) is replaced by a 
Debye screened Coulomb potential, A(q) = g/(q 2 + fi 2 ), and g by the strong coupling 
constant, the above formula leads to da/dt = 47ra 2 /(q 2 + ^i 2 ) 2 . 

One can similarly write down the amplitude for double Pomeron exchange, 

= 2m5(Ek - E f )(-g 2 p ) J 4(P/ - ^ ~ 

where energy conservation at each coupling vertex sets the energy of the internal 
line to £° = y/s/2 = Ef. Amplitudes involving backscattering are suppressed at 
high energies, because of the limited momentum transfer that each the Pomeron can 
impart. If we assume that Pomerons do not overlap in space, the singularity in 
A(q) can be neglected, and the integration over £ z (with respect to the z direction of 
x 2 i = x 2 - xi = Lz + r±) gives us 

M% ] = 2m5(E l -E f )J^u af (p f )(-g 2 p )r (2) u (7i (p i ), 

r (2) = 4(p / -p)X4(p_ p .) e -(p-p i )-x 1 -i(p / -p)-x 2) (6) 



where p = (Ef, yE 2 — £'j_,£±) is the four-momentum of the internal line. One can 
derive the classical Glauber multiple collision cross section from this amplitude by 
averaging and summing over the initial and final state ensemble of the target |27|. In 
the limit of high-energy and small angle scattering, one can neglect the phase factor 
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in the above equation and obtain the amplitude [as defined in Eq. (|2|)], 

/ 2 ( S ,t) = ^/^[-x(b,,)]V^ b , (7) 

where 1/2! comes from the different orderings of the target potentials and = 
P/± — Pi± is the total transverse momentum transfer from the multiple scatterings. 

The eikonal function x CTliCr2 (b, s) is defined as the Fourier transform of the single 
scattering amplitude (besides a factor i/2Ti), 

dtJ (bV^t)T(s,t), (8) 



where t = — and Jo is the zeroth order Bessel function. In the definition of the 
product of eikonal functions, summation over the polarizations of the intermediate 
lines is implied. If we are only interested in unpolarized collisions, we can regard 
T(s, q_i_) as a scalar in the following. One can generalize the double scattering ampli- 
tude to multiple scatterings and sum them together to get the total amplitude, 

00 — i r 1 

/(M) = £/»(M) = ^/ >&£^[-x(b, s )]V<^ b 



n=l 



= — / d 2 b[l - e - x(b ' s) ]e iq ^- b . (9) 

2,71 J 

This is the total elastic quark scattering amplitude in the eikonal formalism. This 



amplitude can also be derived from the scattering theory in quantum mechanics fl39 



where the eikonal function is related to the potential experienced by the scattering 
particle. 

One can consider a pp interaction as multiple quark-quark scatterings and 
derive a similar formula for the total amplitude of an elastic pp collision. The effective 
eikonal function for pp collisions now should be given by 

X (b, 8) = -j=J |^-"^ b t(q ± , s)t(-q ± , s)T(s, q ± ), (10) 

where t(q_i_,s) is the Fourier transform of the quark density distribution inside a 
nucleon and T(s, q_i_) is the amplitude of a quark-quark scattering with one Pomeron 
exchange. Here we neglected the contribution from Reggeon exchanges, so that the 
amplitude T(s,qj_) is purely imaginary p8| (or the eikonal function is real); thus 
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the elastic amplitude of pp collisions, Eq.( |9|), is also purely imaginary. Experimental 
measurements of high-energy pp or pp collisions indeed find the real part of the elastic 
amplitude to be small P0 |. 

With the above elastic amplitude, the differential elastic cross section for pp 
collisions is given by Eq. (01). Using the identity 



/ dtJ (by/^i)J {b'y/^i) = -6(b - b') 
J-oo b 



one can obtain the elastic cross section, 



d 2 b 



1 -e 



-x(M 



(12) 



where we assumed that the imaginary part of the eikonal function is negligible. Using 
the optical theorem, one can also get the total and inelastic cross sections of pp 
collisions | 



Otot 

0"in = O"tot — °cl 



4vrlm/(s,t = 0) = 2 J d 2 b 



1 - e~ x(M) 



d 2 b 



1 _ e -2x(M 



(13) 
(14) 



The approximation of a small imaginary (real) part of the eikonal function 
(elastic amplitude) in Eq. (|TCi| ) also leads us to assume that the quark-quark scattering 
amplitude T(s, q_i_) with one Pomeron exchange is also purely imaginary [cf. Eqs. (Q) 
and fllCP]. According to Eq. (^]) and the optical theorem , T(s, q_i_) should be related 
to the total inclusive cross section of quark-quark scatterings, 



T(s, 0) = -i^-a- mc \(s). 



(15) 



Assuming the dominance of small angle contributions in Eq. (|H]), we have then the 
eikonal function for pp collisions, 



X(b,s) w ^a- mc i(s)T N (b,s), 
where T^(b, s) is the overlap function of two nucleons, 



(16) 



T N {b, s) = J d 2 b't{b', s)t(b - b', s) 



(17) 



which is normalized to one, J d 2 bTj^(b, s) — 1. 
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2.2 Cross Sections: Soft vs. Hard 



To calculate the total inclusive cross section of quark-quark scatterings, 0"i nc i(s), one 
needs to know the Pomeron structure and its coupling to quarks and gluons. In a 
simple model, one can regard a Pomeron as a double gluon exchange with a ladder 
structure as illustrated in Fig. |l|(a). In general, most of the loops inside the ladder 
have soft momenta and thus cannot be calculated via pQCD. We consider this non- 
perturbative contribution to the inclusive cross section as cr so f t (s). As the colliding 
energy y/s increases, contributions from large momentum, pt > Po, to the loop inte- 
gral also increase. As illustrated in Fig. |](b), the cutting diagram of such a Pomeron 
exchange which contains at least one hard loop corresponds to hard scatterings in pp 
collisions. If one assumes that the loop momenta are ordered so that they increase 
toward the hard loop, then the rest of the cut-Pomeron beside the hard loop can be 
resummed to give the quark or gluon (parton) distributions of a nucleon and their 
Gribov-Lipatov-Altarelli-Parisi (GLAP) [fEJ] evolution. The corresponding contribu- 
tion to the inclusive cross section can then be identified as the inclusive jet cross 
section KjJ, 

°jet = L d&dyxdy^: /f jet , (18) 



'2dpj,dy 1 dy 2 



^ K Y, xM*i , P 2 T )x2f b (x2,P 2 T ) dG<X (§ J *' fl) , (19) 



dp\dy x dy 2 ^ dt 

where the summation runs over all parton species, y% and y 2 are the rapidities of 
the scattered partons, and x\ and x 2 are the light-cone momentum fractions carried 
by the initial partons. These variables are related by x\ = xt[& Vi + e V2 )/2, x 2 = 
xr(e~ yi +e~ V2 )/2, xt = 2px/ \/s. The f a (x, Q 2 ) in Eq. fll9|) are the parton distribution 
functions. The pQCD cross sections, da ab , depend on the subprocess variables s = 
Xix 2 s, t = — Pt(1 + exp(y 2 — yi)), and u = — Px(l + exp(yi — y 2 )). A factor K w 2 
is included to correct the lowest order pQCD rates for next to leading order effects 

At small x ~ Pt/Vs, the momentum in the ladder diagram might not be 
ordered anymore and there might be more than one hard loop inside the ladder. 
However, calculations [EH] have shown that such contributions are still small at present 



and future collider energies. If we neglect contributions from multiple hard loops , 
the total inclusive quark-quark cross section can be written as 

Cincl = 0- SO ft(s) + Ojet(s). (20) 
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Hard 



Soft 



(b) Hard = q(x,)q(x 2 ) 



2 




+ q(xi)g(x 2 ) 





?(Xi) g(x 2 ) 



Figure 1: (a) A Pomeron can be divided into the ladder diagrams of a hard part, 
in which at least one loop momentum is large, and a corresponding soft part, (b) 
The hard part of a cut-Pomeron can be factorized into parton distributions and hard 
scattering matrix elements. 



where cx so ft represents the contribution from the cut-Pomeron in which the transverse 
momenta of the loops are limited to px < Po, and therefore is not calculable via 
pQCD. However, as we can see from its definition, cr so ft represents processes involving 
many soft partons in the final state. It cannot be modeled by a simple two parton 
elastic scattering. Since there is no ordering in p?, interference becomes important in 
soft particle production, which can be modeled by, e.g., a string model. 

Given the above decomposition of the total inclusive cross section and the 
eikonal function in Eq. fll6|) , we can rewrite the inelastic cross section for pp collisions, 
Eq. (M), as 



/n. OO 

d 2 b[l - e -(^+^mv(M] = / d%Y,9M ( 21 ) 

J 3=0 

g , {h) = h ct T^(MP e -. jct ^ (M; j > 1; (22 ) 

^(6) = [i_ e -«-«ft^(M] e -«^2V(M)_ (23) 



where gj(b) can be considered as the probability for j hard scatterings [35, 47 1 among 
inelastic collisions at fixed impact paramenter b. Note that a jet in our terminology 
refers to a large pt hard scattering. The average number of jets with px > po 
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in pp collisions is thus (nj e t) = Cjet/^in- The above probabilistic interpretation of 
multiple hard scatterings depends on our assumption that multiple cut-Pomerons 
are independent from each other. This holds as long as the average number of hard 
scatterings is not too large. Given an interaction transverse area, ~ ^/pf], for processes 

? 2 

N ' 



with p T ~ p , independence requires that the total interaction area is less than nRj 



where ~ 0.85 fm is the nucleon radius; i.e., 



crjct ^ (Po-Rv) 2 o" in = a max . (24) 

For p > 2 GeV/c, and cr in ~ 40 mb, the right hand side is cx max ~ 3 barns, and thus 
the independent approximation should hold up to the highest energies foreseen. For 
nuclear collisions the total number of jets is given by 

Ntf = T AA (b)a iet , (25) 

where T AA {b) is the nuclear overlap function at an impact parameter b. For 6 = 0, 
Taa ~ A 2 /ttR\, and multiple mini-jets may be independent as long as 

ojet £ (PoRa) 2 7tR 2 a /A 2 « 2a max /A 2 / 3 . (26) 

For A = 197 the right hand side is 180 mb, and thus independence should apply up to 
LHC energies (see Fig. |2| below). On the other hand, nuclear shadowing of the initial 
structure functions will reduce the jet cross section so that independence should be 
valid beyond LHC energies. 

For p > 1 GeV, o"j et (s) is found to be very small when y/s ^ 20 GeV and 
only the soft component is important. The low-energy data of diffractive nucleon- 
nucleon scatterings exhibit a number of geometrical scaling properties ||48|| in the range 
10 < y/s < 100 GeV, e.g., a e \/a to t = 0.175, and B/a tot = 0.3, where B is the slope of 
the diffractive peak of the differential elastic cross section. This suggests a geometrical 
scaling form |49j for the eikonal function at low energies; i.e., it is only a function 
of £ = b/bo(s), with 7t6q(s) = a so { t (s)/2 providing a measure of the geometrical size 
of the nucleon. We further assume that the nucleon overlap function is given by the 
Fourier transform of a dipole form factor so that, 

T N (b,s) = ^\; (27) 

Xb(0 = §(M) 3 ^ 3 (/ioO, e = &A>(«). (28) 
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Figure 2: The total, inelastic, and elastic cross sections of pp and pp collisions 
calculated by HIJING (solid lines) as compared to the data |58|, |5I], |52|, ^3 
The dashed line is the total inclusive jet cross section with p T > p = 2 GeV/c. 
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where fj,Q = 3.9. With this assumption, the eikonal function can be written as, 
X(b, s) = x(€, s) = Xo(Q[l + ^jct(s)/a soft (s)] . 



(29) 



This form ensures that geometrical scaling^] is recovered at low energies when <Xj Ct «C 

Csoft- 

Choosing p ~ 2 GeV/ c and assuming a constant value of cr so ft = 57 mb at high 
energies, the calculated cross sections and the multiplicity distributions in pp and pp 
collisions agree well with experiments^]. This is the model adopted in HIJING [T^| 
to simulate multiple jet production at the level of nucleon-nucleon collisions. In Fig. 0, 
the calculated total, inelastic and elastic cross sections of pp or pp collisions are shown 



as functions of (solid line) together with experimental data [48], |50|, |5l], [52], H; p4 



The dashed line corresponds to the inclusive jet cross section. The calculated cross 
sections agree well with experiments from ISR to Tevatron and cosmic-ray energies. 
We note that the total inclusive jet cross section increases much faster than the 
inelastic cross section as a function of y/s, leading to an increase of the average 
number of minijets, (Jj et /(r in , with energy. 

We emphasize that the value of p = 2 GeV used in the above calculation is 
only a phenomenological parameter. In order for the model to have some predictive 
power, p should not depend on ^/s. However, its value is subject to considerable 
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controversy |p3|1 . The problem arises from the fact that there is not a clear boundary 
specified by po between soft and hard processes, as shown in Eqs. (|18D and (|20|). If we 
believe that pQCD can be reliably applied to calculate hard contributions to the total 
inclusive cross section up to po, then the rest must be modeled by a phenomenological 
cr so f t . Clearly, cr so f t (s) must depend on the choice of p . With a smaller (larger) p , 
more (less) contributions are included as hard collisions. Hence a so f t (s) would be 
smaller (larger). Obviously many choices of po and cr so f t (s) can give the same total 
cross section cr tot (s). The only restriction is that the sum cr so f t (,s) + o"j et (s) must 
have the right value to give the right energy dependence of the total cross section 
o"tot(s)- Since 0j et (s) increases with decreasing p and <7 so ft(s) is non- negative, po must 
be bounded from below by the experimental data on the total cross section o" tot (s). 
We find that this lower limit in our model is po = 1.2 GeV with the given parton 
distribution functions. For p Q smaller than 1.2 GeV, the inclusive jet cross section 
at high energies is overestimated and the resultant cr to t(s) can never fit the data. In 
addition, one must keep in mind that the lowest po is also bounded by the relevant 
Qo — Po i n the evolution of the parton distribution functions. 

In HIJING, the parton distribution functions are taken to be the Duke-Owens p5| 
parametrization set 1 with Q = 2 GeV. The latest version of this parameterization 
(DOl.l) for f a (x,Q 2 ) is adequate through RHIC energies. However, at higher ener- 
gies, more updated distributions, such as GRV [KB] or MRSD— ' |57f should be used. 



These distributions, constrained by the most recent HERA data joo] in deeply inelas- 
tic ep collisions, are more divergent at small x and give rise to larger minijet cross 
sections, especially at the LHC energy |59[| . Using these distributions, one cannot 
fit (Ttot(s) at high energies with fixed pq. As has been discussed in Ref. J46|, this is 
due to multiple hard loops in one Pomeron exchange. The divergent behavior of the 
distribution functions is usually related to the fact that the hard loops are not ordered 
in pt, which gives rise to multiple jet pair production per cut-Pomeron. Thus the 
inclusive jet cross section in the eikonal function must be normalized by the average 
number of such jet pairs p0|. 



2.3 Modeling the Soft Interactions 

From Fig. |l|(a), we can see that a cut-Pomeron always produces many soft partons 
with small transverse momenta no matter whether there is a hard process present 
or not. Unlike partons from initial and final state radiations associated with a hard 
scattering, these soft partons are not ordered in transverse momentum and coher- 
ence is extremely important, which could virtually produce a coherent color field. 
Therefore, soft interactions may not be modeled simply by regularized parton-parton 
elastic scatterings with small transverse momentum transfers. On the other hand, 
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the color field could also be screened significantly later by the interaction with the 
hard or semihard partons from the hard collisions JB1J in heavy-ion collisions. The 



screening effect will decrease the final particle production from the color field. 

As demonstrated in three-jet events in e + e~ annihilations, the color inter- 



ference effects can be approximated fairly well by a string model |]62 |. In HIJING 
we adopted a variant of the multiple string phenomenological model for the soft in- 
teraction as developed in Refs. |13|, [L4|, Those soft interactions must naturally 



involve small pt transfer to the constituent quarks, as well as induced soft gluon 
radiation which can be modeled by introducing kinks in the strings. The produced 
partons are treated either as hard kinks (for gluons) on the string or the end points 
(for quarks) of another string. The strings are assumed to decay independently via 
quark- ant iquark creation using, in our case, the Lund JETSET7.2 |63|] fragmentation 
routine to describe the hadronization. 

The string excitation is achieved by a collective momentum transfer P = 
(P + ,P~,Pr) between the hadrons. Given initial light-cone momenta 

2 2 
TfL TfL 

Pi = (pf,^r,0 T ), p 2 = (^,p 2 -,O r ), (30) 

Pi P2 

with (pf + rn\/p2){p2 + rn\/pt) = s, the final momenta of the strings are assumed 
to be 



2 2 
TT) TT) 

p[ = (pt-P+^ + P~, PT ), ^ = (^+P + ,ft-P-,- Pr ). (31) 

Pi P2 



The remarkable feature of soft interactions is that low transverse momentum ex- 
change processes with px ^ 1 GeV/c can result in large effective light-cone momentum 
exchanges |13|| , giving rise to two excited strings with large invariant mass. Defining 



P+^^-Hi, p- = x_yS-^, (32) 

P2 Pi 

the excited masses of the two strings will be 

Ml = x-(1-x+)s-Pt, M\ = x+(l - x_)s -p%, (33) 

respectively. In HIJING, we require that the excited string mass must exceed a 
minimum value M cut = 1.5 GeV, and therefore the kinematically allowed region of 
x ± is restricted to 

- x + ) > M* cutl /s, x + (l - x.) > M 2 Tcut2 /s } (34) 
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where M| cutl = M c 2 ut + p 2 Tl M| cut2 = M c 2 ut + p\. Only collisions with 

s > s min = (Mrcuti + M Tcut2 ) 2 (35) 



are allowed to form excited strings. Eq. ( p5|) also determines the maximum that 
the strings can obtain from the soft interactions. In events with both hard and soft 
processes, two strings are still assumed to form but with a kinetic boundary reduced 
by the hard scatterings. 

In HIJING, the probability for light-cone momentum transfer is assumed to 

be 

(1.0 -s ± ) L5 



for nucleons and 



p(x±) = (xi + cys)^ (36) 



P(X±) (x 2 t +C 2 / S ) 1 /4[(l_ x± )2 + c 2/ s ]l/4 ( 3? ) 



for mesons, with c = 0.1 GeV, along the lines of the DPM model |L4|. Soft gluon 
bremsstrahlung processes with p? < Po are also introduced as kinks along the excited 
string. In addition, HIJING also includes an extra low px < Po transfer to the 
constituent quarks and diquarks at the string end points in soft interactions. This 
effect is important at low energies, ~ 20 GeV, to account for the high pj> tails of 



the pion and proton distributions [IS 



I emphasize that the low pt algorithm used in HIJING is a phenomenological 
model needed to incorporate non-perturbative aspects of beam jet physics. Many 
variants of soft dynamics can be envisioned, but none can be rigorously defended 
from fundamental QCD. One of the attractive aspects of going to the highest possible 
collider energies is that the theoretical uncertainties due to soft dynamics are reduced, 
as more and more of the dynamics becomes dominated by calculable semi-hard and 
hard QCD processes. 



2.4 Minijets and Transverse Flow 

To summarize the consistency of HIJING calculations with the available experimental 
data, I show in Fig. ||| the inclusive spectra of charged particles in pp and pp over a 
wide energy range, y/s =50-1800 GeV. We see that the model accounts well for the 
energy dependence of not only the transverse momentum distribution, but also the 
rapidity distribution, as well as the multiparticle fluctuations. This overall quantita- 
tive understanding of multiparticle observables in hadronic interactions, especially the 
magnitude and energy dependence of the conspicuous power-law tail of the pt spec- 
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Figure 3: Data on charged particle inclusive pr spectra Ed 3 a/d 3 p |M[ [65|, |66|. 
pseudorapidity distributions dn/dr] ]57], [58]| , and multiplicity distributions [^, |7D|, |7T| 
P n in pp and collisions compared to HIJING calculations at different energies y/s. 
The transverse momentum spectra and multiplicity distributions have been displaced 
for clarity by extra factors of 10 relative to the absolutely normalized data at y/s = 53 
GeV. 



trum characteristic of pQCD, strongly supports the importance of minijets physics 
at collider energies. To demonstrate the onset of particle production with minijet 
production, I plot in Fig. ^ the energy dependence of the central rapidity density of 
produced charged particles. I also show the contribution from purely soft interactions. 
Since the central rapidity density from soft string fragmentations is almost constant 
as a function of the colliding energy, the increased dn^/dr] mainly comes from the 
hadronization of jets at high energies. The correlation between the central rapidity 
density and minijets is very clear when dn c h/dri is compared with the average number 
of minijets, (rij et ) = Oj Ct /cr in , as functions of */s. 

Beside providing overall agreement with experimental data on multiparticle 
distribution and spectra |T7|, 0, minijet production will also influence correlations 
between observed quantities, some of which could be mistakenly attributed to QGP 
and other collective phenomena. One typical example is the multiplicity and mass 
dependence of the average transverse momentum observed in pp collisions at Tevatron 



Collider energy ||75|| . In hydrodynamic models, transverse collective flow is usually 
generated from the expansion of the thermalized dense system [ 76f . Since all hadrons 
have the same flow velocity, heavy particles tend to have larger transverse momentum 
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Figure 4: dn^/dr]^ = 0) in inelastic pp and collisions as a function of \/s. The 
solid line is the HIJING calculation compared to the data |T7|, |7^, |7^]. The dashed 
line is for events without jet production in HIJING simulations. The dotted line is 
the calculated average number of jets, (rij et ) = Ojet/tfin- 



when they finally freeze out. If the average transverse momentum is plotted against 
the total multiplicity, it is anticipated that for heavy particles it will be larger and 
the increase with the multiplicity will be faster than for light ones. Quite surpris- 
ingly, experiments on pp collisions at Tevatron Collider energy have recently reported 
observation of just this effect J75|] . Levai and Miiller JHJ have studied this reaction 
in a linearized transport theory. They found that there is no time for the baryons to 
equilibrate with the pions during the expansion of a hadronic fireball. They therefore 
suggested that a more novel explanation could be required to account for the appar- 
ent similarity of the flow velocities of the mesons and baryons. They noted that these 
observations could be understood if an equilibrated quark-gluon plasma (QGP) were 



formed in these collisions 78 



However, as also noted by Levai and Miiller, the common transverse flow 
of hadrons may also arise accidentally from the fragmentation of minijets. This is 
indeed what has been found with HIJING calculations [79]. Shown in Fig. |5| is the 



HIJING result of the correlation between (p?) and the total charged multiplicity 
for pions, kaons, and antiprotons (from bottom to top) as solid lines together with 
the data at Tevatron energy. The average multiplicity density (dric^/drj) is calculated 
as rich ( I ^7 1 < A77) divided by 2Ar). The average pr is obtained by applying the same 
procedure as used in the experiment [75[ in which the pt distributions are first fitted 
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Figure 5: (pr) of pions, kaons, and anti-protons (from bottom to top) in —0.36 < 
f] < 1.0 versus the average charged multiplicity density dn^/dr] in \q\ < 3.25. The 
histograms are HIJING results and points are the data ff5|1 . The dashed lines are for 
pp at v/i = 200 GeV. 



with parametrizations [power law a(pT + b)~ n for pions and exponential /3exp(— apr) 
for kaons and antiprotons] and then the fitted parameters are used to calculate {pt) 
in the restricted range < pr < 1.5 GeV/c. Apparently, the data are well accounted 
for by our calculation. Shown as dashed lines in Fig. |5] are the calculated results for pp 
collisions at RHIC energy, ^/s = 200 GeV. They are similar to the results at ^/s = 1.8 
TeV, except that pions have a lower saturated value of (pr) at RHIC energy. Since 
pions are the dominant produced particles, the high multiplicity (p?) f° r ai l charged 
hadrons at y/s = 200 GeV is smaller than at y/s = 1.8 TeV. 

In a model with multiple parton production, it is easy to understand why 
(j>t) increases with n ch . If we decompose the multiplicity distribution into different 
contributions from events with different number of minijets as shown in Fig. |6], we 
find that large multiplicity events are dominated by multiple minijet production while 
low multiplicity events are dominated by those of no jet production. The average 
transverse momentum in events with multiple minijets is certainly larger than those 
without, thus leading to the increase of (pr) with n c h- In order to understand the 
different behavior of the correlations between (pt) and n c h for different particles, we 
recall that the jet fragmentation functions for heavy hadrons tend to be harder than 
for light hadrons, as measured in e + e~ annihilation experiments p0| . Therefore, heavy 
hadrons from jet fragmentation carry larger transverse momenta than light hadrons in 
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Figure 6: Charged multiplicity distributions in NSD pp at y/s = 53 GeV and pp 
collisions at y/s = 200 GeV. The data are from Refs. 7(J . The solid histograms 
are from the HIJING calculation with contributions from events with j = (dot- 
dashed histograms), j — 1 (dashed histograms), and j > 2 (dotted histograms) jet 
production. 



pp or pp collisions. This leads naturally to larger slopes of the (px) vs correlation 
for heavier particles. In other words, the fragmentation of minijets can mimic the 
transverse "flow" effect giving the resultant appearance of collective behavior. 

Now one may ask: if totally different models [77, 81] can both describe the 
data, which one has the true underlying dynamics and how can they be distinguished 
from each other? To answer these questions, one has to turn to the properties of jet 
fragmentation. 



2.5 Resolving Minijets in Hadronic Interactions 

Jets in hadronic interactions, as defined earlier, are produced by large px parton 
scatterings in pQCD. Experimentally, jets are identified as hadronic clusters whose 
transverse energy Et can be reconstructed from a calorimetrical study of the events 
[ |S3|j . However, this cluster-finding method becomes questionable for small and even 
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intermediate Et values due to the background of fluctuations |84| from the underlying 
soft interactions. It has therefore been very difficult to resolve minijets with px ~ 2 
GeV/c from the underlying soft background. 

Although minijets with pr ~ 2 GeV/c are difficult to resolve as distinct 
jets from the background, their effects even in minimum biased events, as we have 
demonstrated, can explain many aspects of hadronic collisions and the associated 
multiparticle production [j33|| . I have proposed that the px dependence of the 
two-particle correlation function can be utilized to study the minijet content in the 
minimum biased events of hadronic interactions. Because particles from jet fragmen- 
tation tend to cluster in phase space, two-particle correlations must be enhanced due 
to minijet production. Especially for two-particle correlations in the azimuthal angle 
<p, contributions from back-to-back minijets should be strongly peaked in both forward 
(A<j) = 0) and backward (A<f) = ir) directions. If we calculate the same correlations, 
but for some selected particles whose transverse momenta are larger than a certain 
Pt cut, the two peaks should be more prominent because these particles are more 
likely to come from minijets. On the other hand, particles from soft production or an 
expanding quark gluon plasma are isotropical in the transverse plane and would only 
have some nominal correlation in the backward direction due to momentum conserva- 
tion. Therefore, the experimental measurement of two-particle correlation functions 
and their p? dependence, especially in pp collisions at the Fermilab Tevatron energy 
where multiplicity and mass dependence of (pt) were first observed, is essential to 
end the present controversy over whether the phenomenon is due to minijets, or string 
interaction [81] or the formation of a quark gluon plasma 77]. 

It is well known that particles from high pt jets are very concentrated in 
both directions of the back-to-back jets. The widths of these high p T jet profiles are 
about 1 in both pseudorapidity rj and azimuthal angle <fi [85]. Minijets, though with 
smaller px, should have similar properties. Since particles with px > PTcut are more 
likely to come from jet fragmentation, we can expect that the two-particle correlation 
functions are more characteristic of jet profiles when pxcnt is larger. This method of 
two-particle correlations is unique because it can determine contributions to particle 
production from minijets which are intangible under the traditional cluster-finding 
algorithm |83| due to their small pt ~ 2 GeV/c. 

The normalized two-particle correlation functions in the azimuthal angle are 
defined as 

/ , , v P(01,02) 1 , Qfi s 

C(01, 02 = TT-s. TT-s ~ 1, 38 

pW>i)pW>2) 

where p(<p) is the averaged particle density in and p(0i,0 2 ) is the two-particle 
density which is proportional to the probability of joint particle production at 4>i an d 
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Figure 7: The correlation functions c(0, A0) vs A0 between two charged particles in 
the full rapidity range in pp collisions at yfs = 1.8 TeV. The dashed line is for all 
charged particles, dotted line for particles with p T > 0.5 GeV/c, dash-dotted line for 
Pt > 1.0 GeV/c, and solid line for p T > 1.5 GeV/c. 



Shown in Fig. |7| are our calculated results on two-particle correlation func- 
tions in pp collisions at y/s = 1.8 TeV for selected particles with different prcut- The 
calculation includes all charged particles in the full rapidity range. As we have ex- 
pected, due to minijets, there is strong two-particle correlation at both A0 = and 
7r, forming a valley at A0 ~ ir/3. For large pxcut> the correlation functions are very 
similar to large px jet profiles as functions of relative to the triggered jet axis flSSfl . 
These features are, however, absent in low-energy hadronic collisions where minijet 
production is negligible |^| . Since there are still many particles from soft production 
which can contribute only to the backward correlation due to momentum conser- 
vation, the study of energy dependence of the relative heights of the two peaks at 
A(f) = and n could provide us information about the energy dependence of minijet 
production. The background at A(f> ~ n/3 also depends on the average number of 
minijets produced|8"2]|. 

Unlike high p? back-to-back jets which are both kinematically bounded to the 
central rapidity region, a pair of minijets can be easily produced with a large rapidity 
gap between them. When we trigger one minijet in a limited rapidity window, the 
other one which is produced in the same parton scattering often falls outside the fixed 
rapidity window. Therefore, if we calculate two-particle correlations for particles in a 
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Figure 8: The same as Fig. [7|, except for charged particles in a limited rapidity range 
of \rj\ < 1. 



7T 



limited rapidity range, the minijet contribution to the backward correlation (A0 
will mostly drop out while the contribution to forward correlation (A0 = 0) still 
remains. Indeed, as shown in Fig. || the forward correlation at A0 = for particles 
in \r)\ < 1 is very strong, but the backward correlation at A<f) = ir is drastically 
reduced as compared to the correlation pattern in the full rapidity range in Fig. [7|. 
Furthermore, due to strong short range two-particle correlations in rapidity ||74|| , the 
forward correlation at A0 = is also enhanced by restricting particles to \rj\ < 1. At 
y/s = 1.8 TeV, we find that the enhancement of backward correlation at A<p = tt due 
to minijets becomes important only when the rapidity window is \rj\ ^ 2. 



3 Parton Production in Nucleus-Nucleus Collisions 

When we extrapolate the model to hadron-nucleus and nucleus-nucleus collisions, 
nuclear effects due to multiple parton scatterings and the interference have to be 
considered. Those nuclear effects can be generally divided into the two categories of 
initial and final state interactions. Initial state interactions and the interference lead 
to an apparent depletion of the effective parton density inside a nucleus, the so-called 
nuclear shadowing. They also lead to a modification of the momentum spectra of 
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produced partons, the Cronin effect |86[, at intermediate energies. However, as I will 
demonstrate below, the Cronin effect will disappear at large transverse momentum 
and at high energies, like other high-twist processes. Both initial and final state scat- 
terings can cause a fast parton to lose its energy. Due to color interference, this energy 
loss is not directly proportional |27|, [31, |32| to the parton's energy. Thus for extremely 
energetic initial partons, this energy loss will be negligible. However, for produced 
partons which travel in the transverse direction, the energy loss becomes important 
relative to their finite transverse momentum. This energy loss will essentially modify 
the parton fragmentation, as well as accelerate parton equilibration. 



3.1 Nuclear Shadowing of Parton Distribution Functions 

Let us start first with the effect of nuclear modification of the parton distributions, 
or nuclear shadowing due to initial state interactions. What we are mostly concerned 
about here is the depletion of the effective parton density which will reduce the initial 
produced parton density and the transverse energy density pi, |87|, H- This Sect 

ion 

is based mainly on the work by Eskola, Qiu and myself in Ref. ||100| . 

"Nuclear shadowing" in the context of the deeply inelastic lepton-nucleus scat- 
tering refers to the measured depletion of the nuclear structure function at small 
Xgj , as compared to Fi of unbound nucleons [jS9fl . The same kind of depletion at 
small x is expected to happen also in the nuclear gluon distributions. During the 
recent years there have been many efforts to explain the measured nuclear shadow- 
ing of quarks and antiquarks |28, 



, [91], [92|, |93], |94| , but for gluons the situation is 
still inconclusive. In these models, shadowing at small x can be attributed to parton 
fusions before the hard scattering which probes the parton distributions. In terms of 
parton fusions, shadowing is also predicted to happen in protons. In this case, "shad- 
owing" refers to the saturation of the actual parton distributions caused by fusions of 
overcrowding gluons at very small x. This mechanism proceeds through perturbative 
Q CD-evolution as formulated in [ 57 , P5| . It has been shown by Collins and Kwiecihski 
that the singular gluon distributions actually saturate due to gluon fusions I will 
first review QCD-evolution of parton distributions and then demonstrate how parton 
shadowing arises from the inclusion of a fusion term in the evolution equation. 

Parton distribution functions inside a nucleon or nucleus are closely related 
to the Pomeron substructure and its coupling to quarks and gluons |37|, [|(J. If one 
assumes that the ladder diagram inside a Pomeron is ordered in px, then one can 
derive a set of evolution equations for the parton distribution functions with respect 
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to the momentum scale of the hard scattering |52|, |97 
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where the splitting functions, P a ^bc(y), are the probability distribution functions 
for the respective radiative processes as illustrated in Fig. |9]( Sbj j x is the fractional 
momentum of the specified partons, and Nj is the number of quark flavors. In 
addition to these splitting processes, one should also take into account the virtual 
corrections as shown in Fig. f|(b). Their contributions to the evolution equations are 
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(41) 
(42) 



These virtual corrections are important to regularize the soft divergences in the split- 
ting functions and guarantee flavor and momentum conservation. If one wants a 
probabilistic interpretation of the evolution and especially a numerical simulation, 
the above virtual corrections will give rise to the Sudakov form factors which are 
necessary to guarantee unitarity. As one can see, the above virtual corrections can 
be incorporated into the evolution equations by introducing a <5-function and '^func- 
tions" (see [^7j for their definitions) in the splitting functions 
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With the above splitting functions one can check that flavor and momentum are 
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Figure 9: (a) Radiative and (b) virtual corrections to parton distributions. 
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In the parton model, a nucleus consists of many valence quarks. Those valence 
quarks, however, are dressed with clouds of sea quarks and gluons due to vacuum fluc- 
tuations (which corresponds to the loops inside a Pomeron ladder). The life time of 
the virtual partons is of the order of the hadron size, Rn ~ 1 fm. If those valence 
quarks are probed by a hard scattering with resolution Q ^> l/Rjy, the coherence of 
the parton clouds will be broken and numerous partons will be released. The parton 
density inside the clouds will then depend on the resolution of the scattering Q. The 
larger the Q the denser the parton clouds, as predicted by the QCD evolution equa- 
tions. These evolution equations are based on perturbative QCD which can only be 
applied down to some scale Qo. Given the initial values of parton distributions at 
Qo, the evolution in Q > Q can be predicted by pQCD. However, the parton distri- 
butions below Qo cannot be obtained from pQCD and so far can only be determined 
from experiments. 
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At small values of x, the leading order QCD evolution equations predict that 
the number of gluons becomes extremely large. It has been known p7| , |95| that for 
sufficiently small values of x and/or of Q 2 , the total transverse area occupied by the 
gluons will be larger than the transverse area of a hadron, so that the interactions 
among gluons can no longer be neglected. Such gluon recombination results in a 
modification of the QCD evolution equations, Eqs. (^) and (0). In the limit of 
small- x and neglecting quark distributions, the modified QCD evolution equation for 



the gluon distribution can be cast in the form |37], |95[ 



d y d t G(y, t) = cG(y, t)- 1 exp(-f - e *) [G(v, t)]\ (49) 

where y = ln(l/x), t = ln[m(Q 2 /A^ CD )], G{y, t) = xg(x, Q 2 ) and c = 12/(11-2^/3) 
with Nf the number of quark flavors. 

The seond term in Eq. ([49]) corresponds to the change of the distribution due 
to gluon recombination. The strength of the gluon recombination is controlled by 
the factor 7, originating from two possible sources. One can consider gluon fusion in 
terms of one Pomeron to two-Pomeron coupling. The two fusing gluon ladders (two 
Pomerons), which couple 4 gluons to 2 gluons, can arise either from independent 
constituents of the proton/nucleus or from the same one, as discussed in []95], |96|, 
I99J. We will refer to the former "independent" and to the latter as "non- 

independent" fusion. Since recombinations from both sources happen simultaneously, 
we divide the parameter 7 into two parts: 

7 = 7i + 7n, (50) 

where 71 corresponds to the independent recombination and 7n to the non-independent 
one. 

To understand the form of the fusion term in Eq. (£49|), one can consider 
xg(x,Q 2 ) as the gluon number n g per unit rapidity (dy = dx/x). If we note that 
the gluon-gluon cross section a ~ a s /Aq CD , the independent gluon fusion probability 
inside a nucleon with transverse area ttR 2 n is then, 

W = ^~-^[xg{x,Q 2 )] 2 - (51) 

This corresponds to the gluon fusion term in Eq. (^9|) (note that e~* ~ a s ). One 
can calculate this fusion term in a QCD recombination model fl95|. In a proton, the 
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strength of the independent fusion then takes the form 



- 2 1 ^ 
7l -3^-2A^ D ' (52) 

where Rn ~ 1 fm is the radius of a proton. 

The non-independent fusion happens inside the same valence quark which is 
assumed to have a transverse size ~ The magnitude of the non-independent 



fusion of the gluon ladders can be estimated as [95 



16 1 vr 3 c 2 
81^27^F'2A^d : 



7ii ~ oT37o77TT2 ' 7772 > ( 53 ) 



where a simplification is made by fixing the initial x of the valence quark to Xi ~ 1. 
We also approximate the scale of the initial valence quark by Qi ~ 2 GeV. 

Let us then consider a large loosely bound nucleus. Naturally, both types of 
fusions are still there but only for the independent one an A^-scaling arises. In this 
case 

A _ 9^4_ tt 3 c 2 

where the nucleus is taken to be a sphere with a sharp surface at Ra = 1.12A 1 / 3 fm. 
The strength of the non-independent fusion remains the same as in the case of a free 
proton: 7^ = 711 . 

It is interesting to notice how the relative contributions of the two types of 
recombination will change when going from a proton to a nucleus of A ~ 200: 711/71 ~ 
7.6 and ^ti/li « 1-0. Thus the non-independent fusion is clearly dominant in a free 
proton while in a large nucleus the contributions from both types are of the same 
order. As a result, parton recombination is strongly enhanced in a heavy nucleus. 

In order to solve Eq. (f49|) exactly by integration, one would need the initial 
distribution either at fixed yo or t and the derivatives along a boundary line (y, t ) or 
(yo,t), respectively. However, since the expression for the non-linear term in Eq. ( pE9|) 
is not valid for the regions where x is large, or where both x and Q are very small, 
the natural boundary condition at x = 1 (or y = 0) is not suitable here. In addition, 
since we do not have sufficient information on other boundary lines, we cannot solve 
Eq. (f49|) by direct integration. Instead, with the semiclassical approximation p7 



we 



use the method of characteristics 10U , so that we can avoid the region of both 
small x and small Q 2 . 

The semiclassical approximation corresponds to neglecting the second order 



29 



derivative term, d y dtln(G), which leaves us with the evolution equation as 

d y z(y,t)d t z(y,t) = c - 7exp[-t - e* + z(y, t)], (55) 

where z(y,t) = \n[G(y, t)]. The reason why this approximation is called semiclassical 
is that it corresponds to using a saddle-point approximation to the integration in the 
integral form of the evolution equation. The above equation can then be cast and 



solved in the form of a set of characteristic equations as shown in detail in |96[ . What 
is needed here is the gluon distribution at the boundary yo[— — ln(x )] at all t, which 
can be obtained by evolving the distribution in the region x > Xq according to the 



original GLAP evolution equations, Eqs. ( |39|) and fl40|) . 

For a proton, we assume that the recombination becomes effective at x ~ Xq ~ 
0.01, which is consistent with |M |99|. We can use the results from a global fitting to 



the parton distributions, like CTEQ ||101|| , to constrain xq. In fact, we will see that 



with Xq = 0.01 the shadowed gluons deviate considerably from the CTEQ gluons only 
after x < 0.001, so our choice for x seems to be reasonable, and we do not expect 
the results to be very sensitive to small changes of x . 

As explained above, the gluon recombination is strongly enhanced in heavy 
nuclei and it starts at somewhat larger values of x than in protons. The corresponding 
boundary line x$ for a nucleus is approximately determined by the relative magnitude 
of the evolution terms in Eq. ([491): for Ga(%o) ~ G( y xo) r fA/l, the relative contribution 



from the gluon fusion in a nucleus is about the same as in a nucleon. This gives x\ 



A 



0.05-0.1. This range of Xq is also supported by other studies [|1 02 1 



Let the total gluonic fractional momentum in the non-shadowed parton distri- 
butions be fo = Jq dx xgcTEq(x, Ql). In the case of a proton, shadowing changes the 
gluonic momentum typically by less than a per cent, which we can clearly neglect as 
a small overall change. 

Perturbative shadowing at small x reduces the gluonic momentum more in a 
nucleus than in a proton. If the momentum fraction of gluons is conserved, there 
must be a corresponding enhancement in the large x region. However, there can 
also be momentum transfers from quarks and antiquarks to the gluons. Here we 
consider nuclei with A ~ 200, for which we expect an overall increase in the glu- 



onic momentum fraction, e^, to be only about 4% [91, |92], 102 1 as constrained by 
the experimental data on quark (anti-quark) shadowing in deep inelastic scatterings. 
These two sources of momentum flow will result in anti-shadowing (enhancement of 
parton densities) at large x > Xq. To account for the anti-shadowing we assume 
g(x,Ql) = a A gcTEQ(x,Ql) for x > Xq with ^cteq(^,Qo) th e non-shadowed gluon 
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distribution and > 1 to be determined by the momentum sum rule, 



f x o f 1 

/ dxxg(x,Q 2 ) + a A dxxg C TEQ(x,Ql) 
Jo c Jx£ 

= fo(l + e A ). (56) 

The right-hand-side of the equation is the gluonic momentum fraction of the non- 
shadowed parton distributions plus the momentum transfer e A from quarks and anti- 
quarks during the recombination. One can solve the above equation iteratively for a A 
with the boundary condition C : g A (x$ ,Qq) = a A g CTE q(xQ > Qo)- Typically, a A ~ 1.1 
for A ~ 200. 

In Fig. |T0|(a), nucleon and effective nuclear gluon distributions for a nucleus 
of A = 200 are compared with the input CTEQ gluon distribution at Q = 2 GeV. 
Notice the ~20 % uncertainty in the nuclear case resulting from varying Xq from 
0.05 to 0.1. To demonstrate the formation of strong perturbative nuclear shadowing, 
corresponding to the relative depletion of gluon distributions in a nucleus, we plot 
the ratio G A (x, Qq)/G(x, Qq) in Fig. |l0|(b). Notice also that as x decreases, the gluon 
distribution in a proton increases much faster, or approaches saturation at a much 
smaller x than that in a nucleus. Therefore, as shown in Fig. |l0|(b) the ratio saturates 
only when the gluons in a proton do so. Thus, saturation of the perturbative nuclear 
shadowing reflects actually the behavior of the gluons in a proton. We see that, due 
to the enhanced gluon recombination in a heavy nucleus, a ~50% nuclear shadowing 
in small- x region is generated perturbatively through the modified QCD evolution, 
accompanied by a ~10 % antishadowing from momentum conservation. 

We have not included the non-perturbative, e.g., Glauber-derived [^8] contri- 
bution to the shadowing. Inclusion of such a contribution is equivalent to changing 
the initial values of gluon distribution at scale Q and will slightly reduce the per- 
turbative contribution to the shadowing through QCD evolution. We also have not 
considered the shadowing of quark distributions in this study of perturbative shadow- 
ing. However, experimental data on quark shadowing are roughly consistent with 
the perturbative gluon shadowing we have just estimated. The shadowing, at least 
for quarks, has been shown to depend weakly on the scale Q \\102 \, also consistent 
with experimental data. 

For a practical implementation of the nuclear shadowing in HIJING, we con- 
sider that quarks and gluons are shadowed by the same amount inside a nucleus and 
use the following parametrized form: 



Af a/N {x) 
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Figure 10: (a) The gluon distributions xg(x, Qq) at Qo = 2 GeV vs. x. The result for 
proton is labeled by xo, and the results for A ~ 200 by x A , respectively. The CTEQ 
gluon distribution ||101|| is labeled by "CTEQ". (b) The ratio xgA^x, Ql)/xg(x, Ql) of 
the shadowed gluon distributions vs. x, demonstrating a strong perturbative nuclear 
shadowing in heavy nuclei. 



1 + 1.19 ln 1/f M [a; 3 - 1.5(>o + x L )x 2 + 3x x L x] 

1.0804 1 / 3 - r 



a A - 



\n(A 

0A(A 1/3 - 1), 



1) 



-x 2 /xl 



(57) 
(58) 



where xq = 0.1 and xl = 0.7. The term proportional to «a in Eq. (pTj) determines 
the shadowing for x < xq with the most important nuclear dependence, while the rest 
gives the overall nuclear effect on the structure function in x > Xq with a very weak 
A dependence. As shown in Fig. , this parametrization reproduces the measured 
overall nuclear effect on the quark structure function in the small and medium x 
regions. However, I should emphasize that this parametrization does not satisfy the 
momentum sum rule and does not include the weak scale dependence which was found 



by Eskola in a detailed study in Ref. ||102 
Eq 

tion 



/] ) represents only the average nuclear dependence of the structure func- 
However, in pA or AA collisions, we have to calculate the effective jet cross 
sections at the nucleon-nucleon level for each impact parameter. Physically, it is nat- 
ural to expect that the nuclear effects on the structure functions could depend on the 
local nuclear thickness at each impact parameter ||87|| . Eq. (|58|) is consistent with the 
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Figure 11: The ratio of quark structure functions Ra{x) = F^ixj/AF^x) as a 
function of x in small and medium x regions for different nuclear mass numbers A. 

and (1581). 



The data are from Ref. [3S] and curves are the parametrization in Eqs. (p 



assumption that the shadowing parameter a^(r) is proportional to the longitudinal 
thickness of the nucleus at impact parameter r. We therefore parameterize the impact 
parameter dependence of a a in Eq. fl5T|) as 

a A (r) = 0.1(A 1/3 - 1)^1 -r'/ift, (59) 

where r is the transverse distance of the interacting nucleon from its nucleus center 
and Ra is the radius of the nucleus. For a sharp-sphere nucleus with thickness function 
Ta{ t ) = 2nFi 2 v 1 — r2 1 R 2 a-> ^ ne averaged a / i( T ') is Jo A Tni^TA^a^r)/^ = a^, with 
a a from Eq. fl58|) . Because the rest of Eq. ([57]) has a weaker A dependence, we only 
consider the impact parameter dependence of a^. 

To simplify the calculation during the Monte Carlo simulations, we can de- 
compose Ra(x, t) into two parts, 

R A (x,r) = R° A (x) - a A (r)R s A (x), (60) 

where a A (r)R s A (x) is the term proportional to a A (T) in Eq. fl5"T|) with aA{r) given 
in Eq. (^) and R A (x) is the rest of Ra^x.t). Both R A (x) and R s A (x) are now 
independent of r. The effective jet production cross section of a binary nucleon- 
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nucleoli interaction in A + B nuclear collisions can therefore be decomposed as 



wn 



o&(r A , r B ) = of* - a A (r A )af ct - a B (r B )af ct + a A (r A )a B (r B )af c ? , (61) 

where crP t , a£ t , a? t and cr^f can be calculated through Eq. fll9f ) by multiplying 
fa(xi,PT)fb(x 2 ,PT) in the integrand with ^(xi)/^^), ^(xi)/?^^), 
-R^(xi)-R^(x 2 ) and -R^(xi)-R^(x 2 ), respectively. 

In central Au + Au collisions at y/s = 200 AGeV, the average parton fractional 
momentum is x = Ip^j^fs ~ 0.02 for mini-jets with > Po = 2 GeV. The impact- 
parameter dependent parton shadowing reduces the averaged inclusive mini-jet cross 
section by 50% from its value in pp. This estimate will be somewhat modified if one 
includes the scale dependence of the shadowing effect. However, at RHIC energies we 
can see from Fig. [TT] that mini-jet production with x ~ 0.02 is still not in the deep- 
shadowed region of x. For sufficiently high energies, most of the mini-jets come from 
x < 0.01 region so that the effective mini-jet cross section may be reduced by a factor 
3 in central Au + Au collisions. Note that this large reduction of mini-jet multiplicity 
by gluon shadowing may increase the limit for independent multi-jet production up 
to \/s = 50 ATeV for Au + Au collisions. 



3.2 Disappearance of Cronin Effect at High Energies 

In hadron-nucleus and nucleus-nucleus collisions, multiple scatterings happen both 
at hadronic and partonic levels. Due to the interference between multiple scattering 
amplitudes, a hadron can only "see" the surface of a nucleus, leading to a total cross 
section ~ R\. Parton fusions which cause nuclear shadowing of the parton distri- 
bution functions can also be viewed as multiple parton scattering in the laboratory 
frame. The parton depletion arises from the destructive interference among different 
scattering amplitudes ||28|| . These multiple scatterings and their interference will also 



affect the momentum spectra of the produced partons. In general, as observed first by 
Cronin et al. [J86| , particle production at large px will be enhanced. However, recent 
experiments |103|1 show that the enhancement decreases both with p? and energy y/s. 



I would like to emphasize that the initial multiple scatterings discussed here 
are very different from the late final parton rescatterings which are responsible for 
parton thermalization. The initial partons in general have very small transverse 



momenta but large center-of-mass (cm.) energies. As demonstrated in Ref. ||27|| , 
multiple scatterings among these initial partons cannot be treated as semiclassical 
cascades which might be valid for final parton rescatterings. The eikonal limit is 
more relevant for the initial partons scatterings which can be treated by a Glauber 
multiple scattering analysis. 
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Following the same procedure that leads to the hadron-hadron scattering cross 
section due to exchange of multiple Pomerons, Eq. (|T4|), one can also obtain the 
inelastic cross section for pA collisions, 

VpA(s) = J d 2 b[l - e -^N(s)T A (b,s)^ ( g 2) 

where <tnn is the total nucleon-nucleon cross section and 

T A (b) = / dzp(b,z), (63) 



is the thickness function of the nucleus at impact parameter b. The nuclear density, 
p(b,z), is normalized to J cPbT^b) = A. 

Let us assume that the total nucleon-nucleon cross section can be described 
by parton scatterings with given cross sections. Define the differential cross section 
for a parton-nucleon scattering, i + N — > j + X, as h^(pi,pj), with the total cross 
section 

^) = ^£^^Wi)- (64) 



Let us also define 



/•±oo 

T^(b,z) = ± / dz'p(b,z'). (65) 

J z 



Note that 



T A {b)=TX(b,z)+T£{b,z). (66) 



Now according to a Glauber multiple scattering interpretation of the pA cross 
section Eq. (|62f) , the probabilities for parton i not to interact with the nucleons up 
to z and for j not to interact after z, are respectively, 

e -<Ti(pi)T-(b,z) ^ e -a-j(pj)r+(M) _ (57) 
Then, the probability for one parton-nucleon scattering, i + A —* j + X, is 



a£1 (i) 
d 2 b 



dze-^ T A^p(b, z)h lj e-^ T i (b ' z) 

oo 

ti j 



(Tj — Oi 

h i3 T A (b) 



-UiT A {b) _ „-vjT A (b) 



l- a -^T A (b) 



+ •■■ , (6* 
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where the definition of the thickness functions T^ib) and Eq. fl66|) have been used, 
and an expansion in terms of aijT^b) has been made in the last step. 

We can use the above result to calculate the probability for a parton k to have 
only one scattering after z, by replacing T^(6) with Tj[(b, z). Then the contribution 
to i + A — > j + X from two parton-nucleon scatterings is 



an {2) 
d 2 b 



/■IAJ _ 

/ dze- a ^^p(b, 



h ik h kj 



-<rkTt(b,z) 



e -°jT+(b,z) 



d 3 p k 
Eh 



E 



d 3 p k h tk h kj 



E k a a 



e -aiT\b) _ e -a k T A (b) 
Ok - &i 



-<TiT A {b) 



e " 3 



UjT A (b) 



(J; 



\n(b)j:^hv+... . 



(69) 



Therefore, up to second order in aTj^(b), we have the cross section for i + A 



J+X, 



^2^t h i kh k j _ {a . + (Tj)h 
. k h k 



i.i 



(70) 



In the second term in the above equation, we can already see that the interference 
between double and single scattering comes into play. If we assume the parton distri- 
bution inside the projectile nucleon to be fi/N^Pi), then the differential cross section 
for N + iV — > j + X is given by 



a 



NN 



E 



d 3 Pi 
Ei 



fi/N(Pi)h ij . 



(71) 



Assuming a hard sphere for the nuclear density p(r) with a radius Ra = -Ro^ 1 ^ 3 ? we 
then obtain the ratio between the differential cross sections of N + A — > j + X and 
N + N -> j + X, 



Ra = 



a 



NA 



Aa NN 



9^/3 i 
16nR 2 ^i 



E 



a NN i 



d 3 Pi 
Ei 



fi/uipt) 



(72) 

Compared to the additive model of hard scatterings, parton production in NA col- 
lisions is enhanced due to multiple scattering. The enhancement is proportional to 
v4 1//3 , which is the average number of nucleons inside a nucleus along the beam direc- 
tion. This formula was used successfully [p.04j] to explain the enhancement of large 
Pt particle production in pA collisions. 

To demonstrate how this enhancement depends on the transverse momentum 
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and colliding energy, let us assume a simple power-law form for the parton-nucleon 
differential cross section, 



h ij = da/dydp 2 T = C/p T , {\y\ < AY/2,p T > Po ). 



(73) 



Assuming that all partons are identical, the total cross section is then, 



1 f , ,2 C tiAYC 
a= 2J dydpT 7 T = (n-2) P g- 2 ' 



(74) 



For large px 3> Po, one can evaluate the integral, 



d 3 p k 
Ev. 



h ik h kj 



dyd 2 qx 



c 



c 



c 



Ft 



-a 



Qt IPt - Qt 

n{n — 2) 



(Po/PtY + Odpo/pr) 



(75) 



Substituting the above cross sections into Eq. ([72]), we have, 

9A 1 / 3 n(n - 2) 



1 + 



a- 



-(Po/ptY + O(( P0 /pt) 



(76) 



There are a few interesting features in the above estimate. The nuclear en- 
hancement of jet cross sections decreases with p T , a general feature of high-twist 
processes. It also depends on the power n of the differential parton cross section. 
This feature is a direct consequence of the interference effect, i.e., the cancellation 
between terms in the single and double scattering amplitudes. From both pQCD 
calculations and experimental data, we know that the power n decreases with the 
colliding energy. Thus, the nuclear enhancement of large p? parton production will 
decrease, and eventually it will disappear at high energies. This trend has already 
been observed in experiments [|103|1 in the energy range y/s = 20 — 40 GeV. A simple 
way to understand this energy dependence is the following. At low energies, the dif- 
ferential cross section for a single hard scattering decreases very fast with P x (large 
power). It is then easier for the incident parton to acquire a large px through two 
successive scatterings, each with small transverse momentum transfer, than through 
a single large px scattering. This is why multiple scatterings in pA collisions cause 
the enhancement of large px parton production. As energy increases, the differential 
cross section for a single scattering decreases less rapidly with p T as compared to low 
energies. It is no longer more economical to produce a large px parton through double 
scatterings than through a single scattering. At extremely high energies, the Cronin 
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effect will disappear. Therefore, we will not consider multiple initial scatterings in the 
following since we only consider parton production in heavy-ion collisions at future 
collider energies (RHIC and LHC). 



3.3 Monte Carlo Simulations 

Since we have argued that multiple initial parton scatterings will become unimportant 
at high energies, we can assume a binary approximation for hard scatterings. In that 
case, multiple hard processes involve only independent pairs of partons. Only very 
rarely does a given parton suffer two high scatterings in one event. As the energy 
increases the number of partons that can participate in moderate pt > Po processes 
increases rapidly and the nuclear shadowing phenomenon will also become important 
at small x. However, for the bulk of parton and transverse energy production, the 
basic independent binary nature of the multiple mini-jet production rate is expected 
to remain a good approximation as long as Eq. (p6|) holds. 

For soft interactions, a nucleus-nucleus collision is also decomposed into binary 
collisions involving in general excited or wounded nucleons. Wounded nucleons are 
assumed to be q — qq string-like configurations that decay on a slow time scale com- 



pared to the overlapping time of the nuclei. In the FRITIOF [13] scheme wounded 



nucleon interactions follow the same excitation law as the original hadrons. In the 



DPM [13] scheme subsequent collisions essentially differ from the first since they are 



assumed to involve sea partons instead of valence ones. In HIJING [|T^] we adopted a 
hybrid scheme, iterating string-string collisions as in FRITIOF but utilizing DPM-like 
distributions as in Eqs. (^) and (|37|). Another difference in the way soft interactions 
are treated in HIJING is that string-string interactions are also allowed to de-excite 
as well as to excite the strings further within the kinematic limits. In contrast, in 
the FRITIOF model multiple interactions are assumed to lead only to excitations of 
strings with greater mass. Many variations of the algorithm for multiple soft interac- 
tion are of course possible as emphasized before. The one implemented in HIJING is 
simply a minimal model which reproduces essential features of moderate energy pA 
and AA data. 

The number of binary collisions at a given nuclear impact parameter is de- 
termined by Glauber geometry. We employ three-parameter Woods-Saxon nuclear 



densities determined by electron scattering data ||105|| to compute that geometry as 



done in ATTILA 15 



For each binary collision, we use the eikonal formalism as given in the previous 
section to determine the collision probability, elastic or inelastic, and the number of 
jets it produces. After a hard scattering, the energy of the scattered partons is 
subtracted from the nucleon and only the remaining energy is used to process the 
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soft string excitation. The excited string system minus the scattered hard partons 
suffers further collisions according to the geometric probabilities. 

For each hard scattering, one then has to take into account the corrections due 
to initial and final state radiations. In an axial gauge and in the leading logarithmic 
approximation, the interference terms of the radiation drop out. The amplitude 
for successive radiations has then a simple ladder structure and the probability for 
multiple emissions becomes the product of each emission |)7| . The virtualities of the 
radiating partons are ordered along the tree, decreasing until a final value /Xq is reached 
below which pQCD is no longer valid. This provides the framework for a Monte Carlo 



simulation of parton showering and its space-time interpretation ||106| , |107 |. 

At a given vertex of the branching tree, the probability for the off-shell parton 
a of virtuality q 2 < g^ ax to branch into partons b and c with fractions z and 1 — z of 



the light-cone momentum is given by ||106| , |107 | 



dq 2 

dV a ^ bc (q 2 ,z) = —dz Pa^bc(z) 



a s [z(l - z)q 2 ]S a {q 2 a ^) 



2tt 



S a { 



where P a _>bc(z) is the Altarelli-Parisi splitting function [42| for the process a 
By requiring the relative transverse momentum qx of h and c to be real, 



(77) 
-> be. 



Qt 



z(l-z 




> 0, 



(7f 



and a minimum virtuality q 2 , q 2 c > *Aq, the kinematically allowed region of phase space 
is then, 



4/io< I 2 <Tn 



e(q)< z <l-e(q), c(q) - -(1 - , ! 



4*18/9*). 



(79) 



Note that the '+function' and 5-function due to virtual corrections in the splitting 
functions in Eqs. (|43|)-(|46|) are not in effect in the allowed phase space. Their role has 
been replaced by the Sudakov form factor S a (q 2 ) which is defined as ||106|, |1(J7 



■-» i f q2 dk 2 /■!-<*) , ^ _ . ,a s \z(l - z)k 2 } . 

S a (q 2 ) = ex V { - I — dzJ2 Pa-+bc(z) l \_ Ji } , (80) 



k 2 Je(k) 



b,c 



27T 



so that n a (g^ ax , q 2 ) = iS a (g^ ax )/5 a (g 2 ) is the probability for parton a not to have any 
branching between g^ ax and q 2 . Therefore, the Sudakov form factor in the Monte 
Carlo simulation is essential to include virtual corrections and ensure unitarity. Since 
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the probability of parton emission between q 2 and q 2 — dq 2 is 



fr\^ Uz) ^m, (81) 

q Je(k) , Z7T 



the probability of no parton emission, by unitarity, will be 

dq 2 /•i-e(fc) .a s [z(l-z)g 21 



d<r /-i-e^j q; s \z 1 — 2)<r , , 



between and g 2 — dg 2 . The probability of no parton emission between g^ ax and 
q 2 — dq 2 then will be the product of the two probabilities, 



n a(gmax> <? 2 ~ d<?) = ^(g^, q 2 ) - dU^q^, q 2 ) 

This is just another form of Eq. fl77|) integrated over z. One can see that the Sudakov 
form factor [|106|| is the solution of the above equation in which unitarity plays an 



important role. 

In principle, one can perform the simulation of initial state radiation processes 
in a similar way. The partons inside a nucleon can initiate a space-like branching 
increasing their virtuality from some initial value Qq. A hard scattering can be 
considered as a probe which can only resolve partons with virtuality up to the scale 
of the hard scattering. Otherwise without the scattering, the off-shell partons are 
only virtual fluctuations inside the hadron and they will reassemble back to the initial 
partons. In PYTHIA, which uses backward evolution, a hard scattering is selected 
first with the known QCD-evolved structure function at that scale, and then the initial 
branching processes are reconstructed down to the initial scale Qq. The evolution 
equations are essentially the same as in final state radiation except that one has to 



convolute with the parton structure functions [[L7| . HIJING explicitly uses subroutines 
from PYTHIA to simulate each hard parton scattering and the associated initial and 
final state radiations. The initial virtuality for the initial state evolution is set to be 
Qo =2 GeV/c, and the minimum virtuality for the final state radiation is /io = 0.5 
GeV/c. The maximum virtuality for the associated radiations in a hard scattering 
with transverse momentum transfer px is chosen to be g max = 2p^. Angular ordering 



can also be enforced in PYTHIA to take into account the soft gluon interference ||107|| 
in the final state radiation. 

After all binary collisions are processed, the scattered gluons from each nucleon 
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are arranged according to their rapidities and connected to the valence quarks and 
diquarks of that nucleon in the collision. The rare hard scatterings of q-q pairs with 
opposite flavors are treated as a special case and processed as independent strings. 

As studied in Refs. [|108| , |109j |110(| , large pr partons must propagate trans- 



versely through the whole excited matter. They will suffer both elastic and radiative 
energy loss. The energy loss dE/dz is very sensitive to the Debye screening mass 
in the medium. Thus, the study of the energy loss of the produced hard jets or jet 
quenching can provide us with a unique probe of the dense matter. To test the sen- 
sitivity of final observables to jet quenching, we used a simple gluon splitting scheme 
(an effective induced radiation) in HIJING, given dE/dz and the mean free path of 
the interaction Xf. The interactions are mostly soft between both soft and hard par- 
tons in the medium. These interactions and induced radiation, in some way, mimic 
pre-equilibrium cascading and semi-thermalization of the produced partons. 

Induced radiative energy loss is modeled in HIJING by determining first the 
points of final state interaction of hard partons in the transverse direction and per- 
forming a collinear gluon splitting at every point. We assume that interactions only 
occur with the locally comoving matter in the transverse direction. Interactions with 
the nuclear fragments are negligible, because the two nuclear discs pass each other 
on a very short time scale IRaIi ^ 1 f m - The interaction points are determined by 
a probability with a constant mean free path Xf, 

dP=?ie- t ' x ', (84) 

where £ is the distance the jet has traveled in the transverse direction after its last 
interaction. 

Since the pre-hadronization state in HIJING is represented by connected groups 
(strings) of valence partons and gluons (kinks), interactions can be easily simulated 
by transferring a part of the parton energy, AE(l) = IdE/dz, from one string configu- 
ration to another. Collinear gluon splitting results in a net jet quenching at the stage 
of hadronization, because the original hard parton energy is shared among several 
independent strings. This simple mechanism of course conserves energy and momen- 
tum and is numerically simple. A more dynamical parton cascade approach involving 
the space-time development of parton showering and collisions between the produced 
partons has been made in PCM |2(| . 



To make sure that the model incorporates the right physics down to low en- 
ergies we have compared the numerical results in pA and AA collisions with data at 
SPS energies. Shown in Fig. [12] are the calculated rapidity distributions of negative 
charged particles in pp (dot-dashed histogram), pAr (dashed histogram) and pXe 
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Figure 12: Rapidity distributions for negative particles in pp (circle, dot-dashed his- 
togram), pAr (crossed-square, dashed histogram) and pXe (square, solid histogram) 
collisions at = 200 GeV. The points are data from Ref. ||111|| and histograms are 
from HIJING calculation. 



(solid histogram) collisions at E\ a ^ = 200 GeV. The data are from Ref. ||1 1 1|| . Because 
jet production is negligible at this energy, particle production occurs mainly through 
soft excitations of projectile and target nucleons. The HIJING low px algorithm re- 
produces an increase of particle production in the central region with the number of 
participating target nucleons. The peak of the rapidity distribution is shifted back 
towards the target region and its height is proportional to the target atomic number. 
In the target region, HIJING under-predicts particle production due to the neglect 
of final state cascading. 

Shown in Figs. ^ and [14] are the calculated rapidity distributions of negative 
particles in central O + Au collisions at = 60 and 200 AGeV, and the transverse 
momentum distributions of negative particles in p + p and central O + Au collisions 
at -Ei a b = 200 AGeV. The overall features of the data ||112| are well accounted for 
except for the enhancement at low px < 0.2 GeV/c in O + Au. That enhancement is 
currently believed to originate also from final state interactions ||113| , |114|| . The data 
for the O + Au collisions are taken with a central trigger. In HIJING simulations, 
central events are selected which can give the corresponding averaged multiplicity. 

The calculated results for central Au + Au collisions at yfs =200 AGeV are 
shown in Fig. [15. The left panel shows the pseudorapidity distributions of charged 
particles. Note that without minijets (dotted), the 2 A soft beam jets in HIJING 
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Figure 13: Rapidity distributions for negative particles in central O + Au collisions 
at E'lab = 60 and 200 AGeV. The data are from Ref. [|112|1 and histograms are from 
HIJING calculation. 
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Figure 14: px distributions for negative particles in pp and central O + Au collisions 
at = 200 AGeV. The data are from Ref. |112| and histograms are from HIJING 
calculation. 
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Figure 15: HIJING results of the dependence of the inclusive charged hadron spectra 
in central Au + Au and p + Au collisions on minijet production (dash-dotted), gluon 
shadowing (dashed), and jet quenching (solid) assuming that gluon shadowing is 
identical to that of quarks and dE/dz =2 GeV/fm with \j =1 fm. R ab (pt) is the 
ratio of the inclusive px spectrum of charged hadrons in A + B collisions to that of 
p + p. 



lead to dN^A/dr] « AdN pp /drj. Each beam jet contributes about 0.75 to the central 
rapidity density almost independent of energy J73J. Without gluon shadowing (dash- 



dotted), minijets are found to approximately triple the rapidity density due to beam 
jets. However, if gluon shadowing is of the same magnitude as that for quarks and 
antiquarks, then the minijet contribution to the rapidity density is reduced by approx- 
imately a half (dashed). The solid histogram shows that the effect of jet quenching 
on the rapidity density is small for dE/dz = 2 GeV/fm and Ay =1 fm. 
Plotted in the middle panel of Fig. ELM is the ratio, 



AB d N AB /dp T /dr] 

R {pt) = <PN„/d£/d v ' (85) 

of the inclusive pt spectrum of charged particles in central Au + Au collisions to 
that of p + p. For > 2 GeV/c, both shadowing and quenching are seen to reduce 
dramatically the inclusive hadron production. In the absence of shadowing and jet 
quenching (dash-dotted) the ratio approaches the number of binary pp collisions in 
the reaction. Shadowing alone (dashed) suppresses moderate px hadrons by a factor 
of about 2. Inclusion of jet quenching with dE/dz = 2 GeV/fm reduces that yield by 
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another factor of about 3-5. It is remarkable and encouraging that the single inclusive 
hadron spectrum appears so sensitive to nuclear effects on jet production. Clearly 
higher order correlations, e.g., high px back-to-back hadron pairs, will be even more 
sensitive to these effects. 

From the above, however, we see that A + A data alone would not be sufficient 
to disentangle the effects of shadowing and jet quenching. To separate the two, 
p + A collisions must be studied at the same energy! In those reactions, the density of 
comoving hadrons is so low that final state interactions leading to jet quenching should 
be negligible. Therefore, any observed suppression of moderate pt hadrons could be 



attributed to gluon shadowing alone. The right panel in Fig. [T5| shows the ratio as 
defined in Eq. (|85|) for central p + Au collisions. Without shadowing (dash-dotted), 
the high px limit again reaches the number of binary collisions, which is about 8 in 
this case. The low pt limit is, on the other hand, controlled by the number of pairs of 
beam jets, which is approximately 4 in central collisions. We see from this figure that 
the ratio R pA (px) at moderate px is indeed sensitive to gluon shadowing. Therefore, 
a systematic measurement of the inclusive hadron spectra in p + p and p + A at the 
same energy is essential in order to determine the magnitude of gluon shadowing. 
Since the Cronin effect on jet production in p + A collisions at y/s ^ f 00 AGeV is 
expected to be small, the spectra of moderate pt hadrons will provide information 
on the gluon structure function complementary to that deduced by the conventional 
direct photon technique that is limited to higher x gluons due to the high tt° —>■ 77 
background. Once the magnitude of gluon shadowing is determined, its contribution 
to the observed suppression of moderate px hadrons in A + A data can be calculated, 
thereby making it possible to isolate the the contribution due to jet quenching. 

3.4 Space-Time Structure of Parton Production 
3.4.1 Space-time history 

So far I have discussed the general features of a pQCD-inspired model for parton and 
particle production in hadronic and nuclear collisions. However, in order to estimate 
the initial number and energy density of the produced partons, one must know their 
formation time. The total formation time must include the interaction time of the 
semihard processes and the formation time for initial and final state radiations. As 
demonstrated in Refs. [EB, E7|, formation time is a consequence of interference in mul- 



tiple scattering processes, which can also be obtained from the uncertainty principle. 
Destructive interference suppresses those secondary scatterings which happen within 
the formation time of the previous interactions. Inclusion of the formation time ef- 
fects is necessary to treat the secondary scatterings properly. Otherwise, one might 
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overestimate the interaction rates and parton production. One source of overestimate 
may arise if one neglects the interaction time of the scatterings. During the inter- 
action time, which depends on whether the scattering is hard, semihard or soft, the 
participating and therefore the produced partons cannot scatter again immediately 
with the beam partons. Another overestimate has to do with the initially radiated 
partons before the hard or semihard scatterings. In calculating jet cross sections and 
the number of hard scatterings, parton structure functions of a nucleus f a /A(x,Q 2 ) 
evaluated at the scale Q 2 = p\ are used. These parton distribution functions then 
have already included the effect of QCD evolution, producing more partons at small 
x and larger Q 2 . Therefore, in Monte Carlo simulations in which initial radiations 
are treated by backward evolution, the radiated partons should not participate in 
any interactions before the end of the corresponding hard scattering. This is also in 
accord with the factorization theorem ||115|| of pQCD. 



Formation time is usually related to a radiative process. According to the 
uncertainty principle, an off-shell parton can be considered as a virtual fluctuation 
and it can only live for a finite time, At, determined by its virtuality q 2 , 

At « q /q 2 , (86) 

where go is the energy of the parton. This At is exactly the formation time of the 
radiation which will follow. Thus, after At, the off-shell parton will then branch or 
"decay" into other partons which can further initiate branchings until a minimum 
virtuality fio is reached. At q 2 < /ig, pQCD is not considered to be valid anymore 
and the process of nonperturbative hadronization takes over. Following this tree of 
branching (which also includes initial space-like radiation) and assuming a classical 
straight line trajectory for partons, one can then follow the space and time evolution 
of the initial parton production ||116 . 



The life time of a virtual parton in Eq. (^) is only a rough estimate according 
to the uncertainty principle. One could also estimate it via a pQCD analysis of the 
radiative amplitudes induced by multiple scatterings |27], |3T|. Thus, qo/q 2 should be 



a good estimate of At in magnitude. However, this formation time or the life time of 
a virtual parton cannot be identified with the decay rate of a virtual parton as has 
been done in PCM [^(J because the later has an additional dependence on a s . The 
decay rate is only related to the radiation intensity, not the formation time of the 
radiated partons. 

Eq. flSTf ) can also be used to estimate the interaction time for each hard scatter- 
ing with q being the sum of the initial or final four-momenta of the colliding partons. 
If the fractional momenta of the partons are X\ and x^, then the interaction time can 
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be estimated as 



In this case, the asymmetric scatterings (x\ 3> X2 or x\ -C X2) have longer interaction 
time than the symmetric ones (xi ~ £2) for fixed parton-parton center-of-mass energy 
X\ X2\fs- We will also assume that the interaction time is the same for all channels. 

In the rest frame of each nucleus, three-parameter Woods-Saxon nuclear den- 
sities are used to construct the nucleon distribution inside the nucleus. The system 
is then boosted to the center-of-mass frame of the two colliding nuclei. Due to the 
fact that gluons, sea quarks and antiquarks are only quantum fluctuations before 
they really suffer scatterings, their longitudinal distribution around the center of the 



nucleon is still governed by the uncertainty principle in any boosted frame ||117|| . One 
usually refers to this distribution as the "contracted distribution" , in which a parton 
with Xi fractional momentum has a finite spatial spread, 

Azi « 2/xisfs. (88) 

Transversely, partons are distributed around their parent nucleons according to the 
Fourier transform of a dipole form factor f35fl . If we define t = as the moment 
when the two nuclei have complete overlap, the overlap of two colliding nucleons can 



happen anywhere inside the diamond-shaped region (with dotted lines) in Fig. [16. 
The diagonal length of the diamond is 2Ra1 with 7 the Lorentz boost factor. The 
interaction point of two partons in a t-z plane can be anywhere within the shaded area 
in Fig. [16|. The solid lines are the trajectories of the two parent nucleons which spread 
around the nuclei according to a longitudinally contracted Woods-Saxon distribution. 

To study the space-time history of parton production, let us look at the results 
of HIJING simulations for Au + Au at the highest RHIC energy, yfs = 200 AGeV. 
Shown in Fig. |Tj] is the total number of produced partons, on-shell as well as off- 
shell, as a function of time. We see that long before the two nuclei overlap and hard 
scatterings take place, partons have already been produced via initial state radiation. 
Note that, if the coherence is not broken by the hard scattering, partons which would 
have been emitted from the initial state radiation will not emerge as produced partons. 



In Fig. [T7| the initiators of the space-like branching are also included as produced 
partons. Therefore, a parton is defined to be a produced one before the hard scattering 
if it has initial state radiation. If a parton does not have initial state radiations, it 



will only become a produced parton after the hard scattering. From Fig. [17], we can 
see that about 2/3 of the total number of partons are produced between t = —0.5 
and 0.5 fm/c while about 200 semihard scatterings happen between t = —0.1 and 0.1 
fm/c as indicated by dashed lines. We find also that about 2/3 of the total number 
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Figure 16: The overlap region in space and time of two incoming partons each with 
spatial spread of Azi and Az2, respectively. The interaction point is chosen randomly 
inside the shaded region. Solid lines show the trajectories of the parent nucleons and 
dotted lines represent the overlap region of the two nuclei. 

of partons are produced in initial and final state radiations. The fraction of partons 
from branching should increase with the colliding energy and with smaller choices for 

fJ,Q. 

To see how hard or semihard scatterings and initial and final state bremsstrahlung 
contribute to parton production, the rapidity distribution of produced partons at dif- 
ferent time is shown in Fig. |l^. Before t — 0, most of the partons come from initial 
state radiation. Since the radiations are almost collinear, these partons move along 
the beam direction and therefore have large rapidities. The semihard scatterings 
then produce partons uniformly over a rapidity plateau and fill up the middle rapid- 
ity region. Final state radiations, which happen after the semihard scatterings, will 
also produce partons uniformly in the central rapidity region. It is quite clear from 
Fig. [18] that the dip of dN/dy in middle rapidity is actually caused by the parton 
production from initial state bremsstrahlung at large rapidity. These partons with 
large longitudinal momenta will move away from the interaction region after the semi- 
hard scatterings. They do not rescatter with the beam partons in the leading twist 
approximation. 

Fig. [19] contains snapshots of dN/ dz at different times to illustrate how the 
parton production evolves in space and time. At t = —0.7 fm/c, as the two nuclei 
approach toward each other before they actually overlap, initial state radiations have 
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Figure 17: The total number of produced partons iV as a function of time t, with 
t = defined as when the two nuclei have complete overlap. 
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Figure 18: The rapidity distribution dN/dy of produced partons at different times t. 
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Figure 19: Parton distribution along the z-axis at different times. 



already begun. These partons have large rapidities and are Lorentz contracted with 
an average spread in z, 



Az ^l/p + 2R A 



2m 



N 



0.25 fm, 



200GeV, 



?9) 



where po = 2 GeV/c is the pt cutoff for semihard scattering, Ra is the nuclear 
radius of Au and is the nucleon mass. After the hard scatterings and during 
the interaction time, partons are produced uniformly in the central rapidity region. 
Afterwards, partons follow straight lines by free-streaming and are distributed evenly 
in z between the two receding pancakes of beam partons (partons from initial state 
radiations). The hyperbolic-shaped equal-time parton distribution in z corresponds 
approximately to a boost invariant parton density in the central region in this free- 
streaming picture. 

Another illustrative way to study the evolution of local parton density is to 
make a contour plot of parton density n in z and t as shown in Fig. ^0[ Here, n is 
defined as 

1 rIN 

(90) 



n 



ttR a dz 



and a sharp sphere distribution is assumed for the nuclear density. One can clearly 
see that partons inside the two approaching nuclei have a spatial spread of Az = 0.25 
fm in z. This spread continues for partons from the initial state radiation as they 
escape from the interaction region along the beam direction with large rapidities. The 
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Figure 20: Contour plot in z — t plane of the parton density n of Eq. fl90"D , as indicated 
by the numbers. The wavy structure along the light-cone is only an artifact of the 
plotting program. 



interaction region where semihard scatterings happen lasts for about 0.5 fm/c, from 
t = —0.25 to 0.25 fm/c. Because of the contracted distribution, some of the partons 
lie outside the light-cone which is defined with respect to the overlapping point of the 
two nuclei. 



If one assumes boost invariance [118 , the parton density can be estimated as 



1 dN 

-kR a t ay 

where r is the proper time and n should be a function of r only. By comparing 
the contour of constant density with the hyperbola of constant r (dot-dashed line) 
in Fig. |20], we see that this is approximately true. We also see that the density, as 
indicated by the numbers, decreases like 1/r due to free-streaming. 



3.4.2 Local isotropy and chemical composition 

Since, as we have demonstrated, there are numerous partons produced within a rather 
short time, the initial parton density is very high at r > 0.25 fm/c, immediately after 
the interaction region (see Fig. ^). Within such a dense system, secondary parton 
scatterings and production are inevitable at a time scale much larger than the average 
formation time. The equilibration time for the system can be estimated by solving a 
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set of rate equations ||119|| as I will show later. In this approach, one must make sure 



that there is approximately local isotropy in momentum space. This could be achieved 
through secondary parton scatterings as has been investigated in Refs. PH|, 120 1. One 



can also use free-streaming to estimate the upper bound of the thermalization time 
[|119| , 121|| , by studying the momentum distribution of partons in a cell of the size 
of the mean free path A/, assuming that parton scatterings are frequent enough to 
achieve thermalization within this upper bound. 

Let's concentrate on the central slice at z = with \z\ < 0.5 fm. At the very 
early stage during the interaction region, produced partons with different rapidities 
are confined to a highly compressed slab with Az pa 0.5 fm. As the system expands, 
partons with large rapidities will escape from the central slice while partons with 
small rapidities remain. Shown in Fig. ^l] is the evolution of momentum distributions 
in p x (solid lines) and p z (dashed lines) at different times (indicated by the number 
on each line). The slope of the p z distribution decreases because partons with large 
longitudinal momenta gradually escape from the central region. At t = 0.7 fm/c, 
the slopes of p x and p z distributions become the same. We should note that the 
isotropy in momentum distribution via free-streaming at this moment is only a tran- 
sient phenomenon which will disappear quickly afterwards. What really determines 
the thermalization time is the relaxation time due to parton scatterings ||122| , |123|| . If 



the relaxation time is short, thermalization can be achieved even before the transient 
momentum isotropy. On the other hand, for a longer relaxation time, the momentum 
distribution will become anisotropic again and will eventually approach thermal equi- 
libration afterwards. To calculate the relaxation time, one has to consider both elastic 
[|122|1 and inelastic |123| parton scatterings and include resummation of hot thermal 



loops. Unfortunately, a complete calculation including higher order processes in a 
nonequilibrium system is not available. In the following we will simply assume that 
the time for local momentum isotropy can be replaced by the time for the central 
slab of about 1 fm to achieve isotropy in momentum space via free-streaming, 

r iso ^0.7 fm/c. (92) 

The readers should be reminded that such assumption will give rise to corrections 
from the neglected viscosities to the results concerning parton equilibration. 

Unlike quarks and gluons in an ideal gas with chemical equilibrium, partons 
produced in the initial stage of heavy-ion collisions, which are determined by the par- 
ton structure functions, the hard scattering cross sections and the radiation processes 
in pQCD, are far from chemical equilibrium. Due to the difference in the numbers 
of degrees of freedom in color space, cross sections involving gluons are always larger 
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Figure 21: p x and p 2 distributions at different times for partons in the central slice 
\z\ < 0.5 fm. 



than those of quarks. Similarly, both initial and final state radiations produce more 
gluons than quarks and anti-quarks. Therefore in pQCD, the ratio of produced quarks 
to gluons is much smaller than the ratio in an ideal gas, which is 9/4 for three quark 
flavors. 



Shown in Fig. [22| are the fractions of produced quarks and anti-quarks as 
functions of time t in pp and AA collisions at y 7 ? = 200 AGeV. For Au + Au collisions, 
hard scatterings with p T > 2 GeV/c produce about 13% quarks and anti-quarks. If 
initial and final state radiations are not included the ratio jumps to 18%, because 
radiations produce more gluons than quarks and antiquarks. In pp collisions, about 
28% of the partons produced via hard scatterings without radiations are quarks and 
anti-quarks. 

The difference between Au + Au and pp collisions is due to the different A 
dependence of the valence and sea quark production. Since the partons are produced 
via binary collisions, the number of produced gluons and sea quarks and antiquarks 
scales like A 4 / 3 . On the other hand, baryon number conservation requires valence 
quark production to scale like A. Given the fraction of total produced quarks q pp = 
0.28 and valence quarks v pp = 0.14 in pp collisions, one can find for Au + Au collisions, 



Qaa 



Qpp - (i - A-vzy 

1 - (1 - A~V 3 )v. 



(93) 
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Figure 22: The fractional number of produced quarks and anti-quarks as a function 
of time for Au + Au (solid) collisions, Au + Au without radiations (dashed), and p + p 
without radiation (dot-dashed). 



which gives q^A — 0.18 for A = 197, as we obtained from Fig. |22| of the numerical 
calculation. To demonstrate the effective valence quark production, we plot in Fig. ^ 
the rapidity dependence of the fractional quark number. In pp collisions, the valence 
quark production peaks at large rapidities. In nucleus-nucleus collisions, the valence 
quark production has a different scaling in A than gluons and sea quarks. This is why 
the fraction of quark and antiquark production is suppressed more at large rapidity 
than in the central rapidity region. The relative quark and antiquark production in 
the central rapidity region is further suppressed by final state radiation. 



The relative quark production found in PCM |124fl is about five times larger 
than what I show here. This is due to an overestimate of intrinsic quark production 
through flavor excitations in PCM. Such flavor excitations should be suppressed by 



interference among pQCD amplitudes at the same order | 125 |. For the same reason 



other heavy quark production is also overestimated in PCM ||126|| . 

The small number of produced quarks relative to gluons within pQCD has 
important consequences for chemical equilibration of the partonic system. Because 
of the small initial relative quark density and small quark production cross section 
as compared to gluon production, it takes a very long time for the system to achieve 
chemical equilibrium, if it ever does. If this time is longer than the phase transition 
time, a fully equilibrated QGP may never be formed [ |119|| . 
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Figure 23: The final fractional number of produced quarks and anti-quarks as func- 
tions of rapidity for Au + Au (solid) collisions, Au + Au without radiations (dashed), 
and p + p without radiation (dot-dashed). 



3.4.3 Multiple initial parton scatterings ? 

I have already demonstrated earlier in Sec. 3.2 that multiple parton scatterings be- 
come less important at high energies and eventually will disappear. This phenomenon 
can also be understood in a space-time picture. But before doing that, I should em- 
phasize again that the multiple parton scatterings in this context always refer to those 
during the overlapping time of two colliding nuclei. They are different from the final 
parton rescatterings which are responsible for parton thermalization. 

We have explicitly taken into account the interaction time for the semi-hard 
scatterings, which is roughly Tj ~ 1/pr- Inside a highly Lorentz contracted nucleus, 
the spatial spread for partons which could participate in a hard scattering with trans- 
verse momentum transfer p? is also about Az ~ 1/pr- This would leave the produced 
partons no time to have another hard scattering of pt with the incoming beam par- 
tons. For finitely contracted nuclei at relatively low energies, this kind of double 
high pj- scattering is still possible. However, these higher twist processes should be 
suppressed by a factor of \ jp\- It is also possible for a parton to go through a hard 
and a soft scattering subsequently since the soft partons always have a spatial spread 
of 1 fm. This kind of hard-soft multiple interactions constitutes the leading contribu- 
tion to higher twist corrections to hard processes in nuclear collisions ||127|| . However, 
double semihard scatterings at high energies with p T ~ p will be suppressed due to 
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finite interaction time. 

The large pt enhancement of both Drell-Yan dilepton production and single 
hadrons in pA collisions at y'i < 50 GeV is considered as a result of multiple parton 
interactions |[L28|1 . At these energies, the contracted length of a heavy nucleus is 



still relatively larger than the interaction and formation time. The factorized form 
of the parton model is then modified due to the finite beam energy. However, a 
collection of experimental data |103|1 even in this intermediate energy range has 



already shown the effects of the finite interaction time. As the energy increases, the 
interaction time becomes more important as compared to the size of an increasingly 
contracted nucleus. The partons then have less time for secondary scatterings. This 
then leads to the observed decrease of the large px enhancement in pA interactions as 
the beam energy increases. At ultrarelativistic energies, one should therefore expect 
the suppression of multiple semihard scatterings. As I have argued earlier, partons 
from the associated initial state radiation will not scatter again with the incoming 
beam partons due to factorization. 

4 Parton Equilibration 

The initial semihard interactions between two beams of partons happen in a short 
period of time. After the leading partons leave the interaction region, a partonic 
system will be formed in the central region. In this section, I will discuss how this 
partonic system equilibrates and describe the physics involved, based on a microscop- 
ical description. The thermalization and equilibration in this description are driven 
by multiple scatterings and parton production (or annihilation) by induced radiation 
(or absorption). Though an exact quantum field treatment of multiple parton inter- 
action is not possible with present technology, a semiclassical approximation has been 
undertaken ||20|| . The most important and difficult effort in a semiclassical approxi- 



mation is to simulate the quantum effects, like interference, coherence, and especially, 
nonperturbative phenomena in QCD. Under extreme conditions at high temperature 
and densities, these quantum effects become very important and may dictate the evo- 
lution of the partonic system. In this section, I will also discuss the interference effects 
associated with induced radiations, since radiative processes produce extra partons, 
leading to fast parton equilibration. 
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4.1 Initial Conditions: a Hot and Undersaturated Gluonic 
Gas 



As we have discussed earlier, the kinematic separation of partons with different rapidi- 
ties in a central slab of about 1 fm give rise to a transient local momentum isotropy 
at the time of the order of r iso = 0.7 fm/c. We will further assume this is the ac- 
tual kinetic equilibration (or thermalization) time for the partonic system. Though 
this assumption is not well founded, this is the best we can do before a complete 
calculation of the relaxation time due to parton rescatterings and radiations is avail- 
able. After this time the momentum distribution of partons is assumed to be locally 
isotropic and approximately exponential. What we are interested next is how this 
system approaches chemical equilibrium. 

Let us first estimate the initial conditions for the equilibration from the previ- 
ous chapters. Since we are here primarily interested in the chemical equilibration of 
the parton gas, we shall assume that the parton distributions can be approximated 
by thermal phase space distributions with non-equilibrium fugacities A^: 



f(k;T,X l ) = \ l (e l3u - k ±\ l ) , (94) 

where (3 is the inverse temperature and is the four-velocity of the local comoving 
reference frame. When the parton fugacities Aj are much less than unity as may 
happen during the early evolution of the parton system, we can neglect the quantum 
statistics in Eq. fl94|) and write the momentum distributions in the factorized form, 

f(k;T,X i )=X i (e^ k ±i)~\ (95) 

Using this form for the distributions, one has the parton and energy densities, 

n = (A 9 ai + Ag&OT 3 , s = (X g a 2 + X q b 2 )T 4 . (96) 

where a± = 16C(3)/7r 2 « 1.95, a 2 = 87r 2 /15 ~ 5.26, for a Bose distribution, b\ = 
9C(3)iV / /7r 2 » 2.20 and b 2 = 77V 2 N f /20 » 6.9 for a Dirac distribution for a baryon 
symmetric system, X q = Xq. Since we have demonstrated that boost invariance is 
a good approximation for the initially produced partons, we can then estimate the 
initial parton fugacities, A° and temperature T from [cf. Eq. (|9~1|)1 



1 dN 4 
7rR 2 A T iso dy ' £ ° H ° n 



no = —^2z — ~ ' e = n -(k T ), (97) 



where [kr) is the average transverse momentum. The quark fugacity is taken as 
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RHIC LHC 

r iso (fm/c) 07 0.5 
e (GeV/fm 3 ) 3.2 40 
n (far 3 ) 2.15 18 
(k T ) (GeV) 1.17 1.76 
T (GeV) 0.57 0.83 
A° 0.05 0.124 
Aj 0.008 0.02 



Table 1: Values of the relevant parameters characterizing the parton plasma at the 
moment r iso , when local isotropy of the momentum distribution is first reached. 



A° = 0.16A°, corresponding to a ratio 0.14 of the initial quark(antiquark) number to 
the total number of partons. Table [I] shows these relevant quantities at the moment 
Ti SO , for Au + Au collisions at RHIC and LHC energies. One can observe that the 
initial parton gas is rather hot, reflecting the large average transverse momentum. 
However, the parton gas is very far from saturation in phase space as compared to the 
ideal gas at the same temperature. The gas is also far away from chemical equilibrium 
with respect to parton species (gluons and quarks with different flavors) since it is 
dominated by gluons. I should emphasize that the initial conditions listed here result 
from HIJING calculation of parton production through semihard scatterings. Soft 



partons, e.g., due to parton production from the color field |nj], are not included. 



4.2 Master Rate Equations 

In general, chemical reactions among partons can be quite complicated because of 
the possibility of initial and final-state gluon radiations. As we have discussed earlier, 
interference effects due to multiple scatterings inside a dense medium lead to a strong 
suppression of soft gluon radiation. One lesson we learned from the LPM effect is 
that the radiation between two successive scatterings is the sum, on the amplitude 
level, of both the initial state radiation from the first scattering and the final state 
radiation from the second one. Since the off-shell parton is space-like in the first 
amplitude and time-like in the second, the picture of a time-like parton propagating 



inside a medium in the parton cascades simulations |20| breaks down. Instead, we 
shall here only consider processes where a single additional gluon is radiated, such 
as gg — ► ggg between two scatterings, in which we can include the LPM effect by 
a radiation suppression factor. The analysis of QCD LPM effect in Refs. |27|, |31| 
has been done for a fast parton traveling inside a parton plasma. We will use the 
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results here for radiation off thermal partons whose average energy is about T, since 
we expect the same physics to happen. 

In order to permit the approach to chemical equilibrium, the reverse process, 
i. e., gluon absorption, has to be included as well. This is easily achieved by making use 
of detailed balance. Closer inspection shows that gluon radiation is dominated by the 
process gg — > ggg, because radiative processes involving quarks have substantially 
smaller cross sections in pQCD, and quarks are considerably less abundant than 
gluons in the initial phase of the chemical evolution of the parton gas. Here we 
are interested in understanding the basic mechanisms underlying the formation of a 
chemically equilibrated quark-gluon plasma, and the essential time-scales. We hence 
restrict our considerations to the dominant reaction mechanisms for the equilibration 



of each parton flavor. These are just four processes 129 



99 +-> 999, 99 <-> 99- (98) 

Other scattering processes ensure the maintenance of thermal equilibrium (gg «-» 
99-, 91 9Qj e ^ c or yield corrections to the dominant reaction rates (gq <-> qgg, 
etc.). 

Restricting ourselves to the reactions Eq. (|98"D, the evolution of the parton 
densities is governed by the master equations: 

= ^(1-^1 (99) 




^(v") = ^K^)^(i-^0], (ioo) 

where hi = rij/Aj denote the densities with unit fugacities, Aj = 1, <7 3 and <r 2 are 
thermally averaged, velocity weighted cross sections, 

^3 = (o{gg -> ggg)v), o 2 = (a(gg qq)v). (101) 

We have also assumed detailed balance and a baryon symmetric matter, n q = rig. 
If we assume that parton scatterings are sufficiently rapid to maintain local thermal 



equilibrium, and therefore we can neglect effects of viscosity due to elastic [61], [L30 



and inelastic ||123|| scatterings, we then have the hydrodynamic equation 



^(su") + P OX = 0, (102) 

which determines the evolution of the energy density. For a more complete study, one 
can include the viscosity corrections to the energy-momentum tensor and modify the 
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above hydrodynamic equation. This is beyond the scope of this review. 

For a time scale r -C Ra, we can neglect the transverse expansion and consider 
a purely longitudinal expansion of the parton plasma, which leads to Bjorken's scaling 
solution [|118|| of the hydrodynamic equation: 



£ + = 0. (103) 

GST T 

We further assume an ultrarelativistic equation of state, e = 3P, with rij and 
e given by Eq. (|96p. We can then solve the hydrodynamic equation, 

[A s + — AJ 3/4 T 3 r = const. , (104) 

and rewrite the rate equation in terms of the time evolution of the parameters T(t), 
X g (r) and A 9 (r), 

£+*f + 7 = fl3 (l-A 9 )-2fl 2 (l-|) (105, 




^ + 3^ + 1 = r£L[?-?), (106) 

where the density weighted reaction rates i?3 and R2 are defined as 

#3 = l^n g , R 2 = \02ng. (107) 

Notice that for a fully equilibrated system (A fl = \ = 1), Eq. ( |104| ) corresponds to 
the Bjorken solution, T(r)/T = (to/t) 1 / 3 . 

4.3 Parton Equilibration Rates 

Often referred to as Landau- Pomeranchuk-Migdal (LPM) effect, the interference in 
radiation induced by multiple scatterings was first investigated in QED for induced 
photon radiation off a charged particle propagating through a dense medium [^9], |30 



It suppresses those soft photons whose formation time is much larger than the mean 
free path of the fast particle inside the medium. In QCD, however, the interference 
pattern is different from the QED case, since multiple scatterings through gluon 
exchanges will force a quark to change its color along its path. In recent studies 
27], [TT], pQCD analyses indeed indicated that the LPM effect suppresses soft 



gluon radiation by multiple scatterings. Due to color exchanges, the interference has 
some non-Abelian behavior which actually depends on the color representation of the 
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jet parton. The same mechanism also causes the color conductivity in a quark-gluon 



plasma to be small [131, 132 



To account for the LPM effect in the calculation of the reaction rate R3 for 
99 999-1 we w iH use the triple differential cross section from Refs. [27|, ^TJ and apply 



the LPM suppression to the radiation of gluons whose effective formation time tqcd 
are much longer than the mean free path A/ of multiple scatterings. At the same 
time, the LPM effect also regularizes the infrared divergency associated with QCD 
radiation amplitude. However, er 3 still contains infrared singularities in the gluon 
propagators associated with QCD scatterings. For an equilibrium system one can in 



principle apply the resummation technique developed by Braaten and Pisarski ||133 
to regularize the electric part of the propagators, though the magnetic sector still has 
to be regulated by an unknown magnetic screening mass which can only be calculated 



nonperturbatively ||134|| up to now. Since we are dealing with a nonequilibrium sys- 



tem, Braaten and Pisarski's resummation is not well defined. As an approximation, 



we will use [135] the Debye screening mass, 



= 4" r k f( k ) dk = 4rra s T 2 \ g , (108) 
ix Jo 



to regularize all singularities in both the scattering cross sections and the radiation 
amplitude. 

To further simplify the calculation we approximate the LPM suppression factor 



|HJ by a ^-function, 6(Xj — tqcd)- The modified differential cross section for gg — > ggg 
da* da 9 J dn n 



is then 

' 7? - <UT ' A ""* -^-TQcoMv^-^coshy), (109) 



dq]_dyd 2 kx dq\ dyd 2 k\ 



where 



C A 2 cosh y 

TQCD = 2c 2 ^T~ (110) 



is the effective formation time [27, |3T|]. The second step- function accounts for energy 
conservation, and s = 18T 2 is the average squared center-of-mass energy of the two 
gluons in the thermal gas. The regularized gluon density distribution induced by a 
single scattering is, 

dn g _ C A a s ql 



dyd 2 k ± vr 2 fel[(kj.-qj.) 2 + ^]" 
Similarly, the regularized small angle gg scattering cross section is 



doff _ 9 2na 2 s 

2\2' \ llZ ) 



dq 2 ± 4 (ql + fi 2 D ) 
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The mean free path for elastic scatterings is then 



rs/4 fid" 9 1 

*f = n*L d ^ = 8 aiasT l + 8na s \ 9 /9 ' (U3) 



o 



which depends very weakly on the gluon fugacity X g , where again at = 16C(3)/7r 2 
1.95. Using 



27r . , 1 2tt 



(k ± -qu) 2 +^ J(kl + q 2 ± +fj, 2 D ) 2 -Aqlk 2 ± 



(114) 



we can complete part of the integrations and have 

32 

Rs/T = —a s X g (l + 87ra s A 5 /9) 2 X(A 9 ). (115) 
6<X\ 

Here, X(A S ) is a function of A g : 

, rV**f , rs/^n , z \ cosh _1 (yi) 

I(X g ) = dx dz j v 



(1 + ^) 2 \x^[x + (I + z)x D } 2 - Ax z 



XD 



+ 4t ; cosh " 1 (v /i ) ^ (us) 

sA / ^[1 + X (l + 2;)y ]2 - 4x 2 ?/ D 

where xd = I^dX 2 and yo = l^bl s - We can evaluate the integration numerically and 
find out the dependence of R3/T on the gluon fugacity X g . In Fig. [24], R3/T is plotted 
versus A ff for a coupling constant a s = 0.3. The gluon production rate increases with 
X g and then saturates when the system is in equilibrium. 

The calculation of the quark equilibration rate R2 for gg — > qq is more straight- 
forward. But one also has to use effective quark and gluon propagators to regularize 
the singularities for massless quark production [|136j| . We can also estimate the quark 



equilibration rate by using the effective thermal quark mass [137 



( A„ + -A (/ ) ^-a s T 2 , (117) 



as a cutoff in the divergent integral over momentum transfer. After averaging over 
the thermal gluon distribution, we have 119| 

R2 = \e2n g w 0.24iV) a s %Tln(5.5/A g ). (118) 
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Figure 24: The scaled gluon production rate R3/T (solid line) for gg — > ggg and the 
quark production rate R2/T (dashed line) for gg — > qq are shown as function of the 
gluon fugacity X g for ot s = 0.3. 

The dashed line in Fig. shows the normalized rate R2/T for Nf = 2.5, which takes 
into account the reduced phase space of strange quarks at moderate temperatures, as 
a function of the gluon fugacity. 

4.4 Evolution of the Parton Plasma 

With the parton equilibration rates which in turn depend on the parton fugacity, we 
can solve the master equations self-consistently and obtain the time evolution of the 
temperature and the fugacities. Shown in Figs. |25] and [26| are the time dependence of 
T, X g , and X q for the initial conditions listed in Table [l] at RHIC and LHC energies. 
We find that the parton gas cools considerably faster than predicted by Bjorken's 
scaling solution (T 3 r = const.) shown by the dotted lines, because the production 
of additional partons during the approach to chemical equilibrium consumes an ap- 
preciable amount of energy. The increased cooling, in turn, impedes the chemical 
equilibration process; this is more apparent at the RHIC than at the LHC energies. 
Therefore, the parton system can hardly reach its equilibrium state before the effec- 
tive temperature drops below T c « 200 MeV in the short period of time of 1-2 fm/c at 
the RHIC energy. At the LHC energy, however, the parton gas, especially the gluon 
component, becomes very close to equilibrium (in terms of phase space occupation) 
and the plasma may exist in a deconfined phase for as long as 4-5 fm/c. Another 
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Figure 25: Time evolution of the temperature T and the fugacities X g and X q of gluons 
and quarks in the parton plasma created in Au + Au collisions at the RHIC energy 
of -y/s = 200 AGeV. The initial values for T. 
simulations and are listed in Table [TJ. 



X g and X q are determined from HIJING 



important observation is that quarks never reach chemical equilibrium at both ener- 
gies. This is partly due to the small initial quark fugacity and partly due to the small 
quark equilibration rate. 

I should emphasize that the initial conditions used here result from the HI- 
JING model calculation in which only initial direct parton scatterings are taken into 
account. Due to the fact that HIJING is a QCD motivated phenomenological model, 
there are some uncertainties related to the initial parton production. 

(1) As I have discussed in the second section, the total jet cross section and 
initial parton production depend on the pr cutoff po- We only included those partons 
which are produced through a hard or semihard scattering with transverse momentum 
transfer px larger than pq. Soft partons through soft interactions with px < Po are 
locked up in the phenomenological strings which will fragment into hadrons later. 
How one should include these soft partons in the parton equilibration processes and 
what value one should choose for po cannot be answered within the model. The initial 
parton density is especially sensitive to the value of Pq. 

(2) The dominant initial gluon production is also sensitive to the gluon distri- 
bution at small x. But so far there is no precise measurement of the gluon distribution 
in the small x region where most of the minijets originate at RHIC and especially 



at the LHC energies. The constraints from recent HERA data [58] suggest a more 
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Figure 26: The same as in Fig. [25|, except for LHC energy, ^fs = 5.5 ATeV. 



singular behavior of parton distributions than the original Duke-Owen parametriza- 
tion which we have used here. A study in Ref. |59 shows that these more singular 
gluon distributions can give rise to a much larger amount of produced partons at the 
LHC energy. Furthermore, we have no handle experimentally so far on the nuclear 
shadowing of the gluon structure function. In HIJING we have assumed the same 
shadowing for gluons and quarks, which has reduced the parton production by about 
half inside a heavy nucleus. 

(3) Secondary elastic parton scatterings before r iso have not been considered 
here, which could speed up the thermalization process. For example, the parton 
cascade model f2(| claims a shorter thermalization time than what we have used here. 
Changes in the initial thermalization time, therefore, will also lead to uncertainties 
in the initial energy and parton densities. 

We can estimate the effect of the uncertainties in the initial conditions on the 
parton gas evolution by multiplying the initial energy and parton number densities 
at RHIC energy by a factor of 4. This will result in the initial fugacities, A° = 0.2 



0.024. With these high initial densities, the parton plasma can evolve 



and A° 

into a nearly equilibrated gluon gas, as shown in Fig. ^27\. The deconfmed phase 
will also last longer, for about 4 fm/c. However, the system is still dominated by 
gluons, and thus has fewer quarks and antiquarks than expected in a fully chemical 
equilibrated system. If the uncertainties in the initial conditions are caused by the 
soft parton production from the color mean fields, the initial effective temperature 
will decrease because soft partons are generally less energetic. Therefore, we can 
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Figure 27: The same as in Fig. except that the initial parton densities are 4 
times higher than given in Table |T] with the same (ordinary lines), or reduced initial 
temperature, T = 0.4 GeV (lines with stars) 



alternatively increase the initial parton density by a factor of 4 and decrease T to 
0.4 GeV at the same time. This leads to higher initial fugacities, X° g = 0.52 and 
\q = 0.083. As shown in Fig. |27| by the curves with stars, this system evolves faster 
toward equilibrium, however, with shorter life-time in the deconfmed phase due to 
the reduced initial temperature. We thus can conclude that perturbative parton 
production and scatterings are very likely to produce a quark-gluon plasma (or gluon 
plasma to be more accurate) in ultrarelativistic heavy-ion collisions at LHC energy. 
However, the fate of the quark-gluon plasma at RHIC energy has to be determined 
by a more careful examination of the uncertainties in the initial conditions. 

4.5 Effects of gg —> [n — 2)g processes 

So far we have only considered gg — *> ggg process in parton equilibration. In fact, 
higher order processes like gg — > (n — 2)g for n > 5 should also contribute, where n is 
the total number of gluons involved in the processes. Even though these processes are 
higher order in a s , the infrared behavior of QCD radiation gives rise to a logarithmic 
factor In Q 2 after integration which can offset the logarithmically decreasing coupling 
constant and make the contributions relatively important for all orders in a s (Q). This 
is especially true for interactions with large momentum transfer Q. 

The complete amplitude involving n-gluon processes in general has contribu- 
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tions from different helicity transitions. Parke and Taylor ||138|| have derived the exact 
form of the maximum helicity violation amplitude for the n-gluon processes, 



— 1 f>J . p S12S23 ■ ■ ■ S-nl 



where Sij = {pi+Pj) 2 and the summation P is over the (n — 2)!/2 non-cyclic permuta- 
tion of (1 • • -n). There is no known exact expression for the other chiral amplitudes. 
Kunszt and Stirling |139|1 proposed to add a factor to the Parke-Taylor amplitude, 

\M* S \ 2 = C KS (n)\M p J\ 2 , with C KS (n) = 2 "~ 2(n i H : 1) , (120) 

n[n — 1) 

to get the total n-gluon amplitude. This procedure can be checked to recover the 
analytic results for n = 4 and 5 ||140| . 



Using the above tree-level n-gluon amplitude, Xiong and Shuryak ||141|| have 
calculated the n-gluon reaction rate, 

o 1 f d 3 h d 3 k 2 . 

K gg ^(n-2 )g - 2 , (n _2)\J (27r)3 2c 1 /lfclj (2vr)3 2^ /(fc2j 

\MT\\^f5\k X +k 2 -Y: ki) n[i + /(^j^A-. (121) 

i=3 i=3 \ / 1 

In order to regulate the infrared and collinear divergencies they impose a cut-off, 

>s = r? 2 T 2 , (122) 

in the integration over the phase space. If such a cut-off is provided by the Debye 
screening mass, then rf = a s /4n. Neglecting the temperature dependence of the 
coupling constant, the rates scale with temperature as 

R gg ^(n-2) g = a n (?7)T 4 , (123) 

where the dimensionless coefficients a n (n)'s depend only on the screening parameter 
rj. Xiong and Shuryak found that for 77 = 1.1, oi n (rj) decreases with n. They also 
found that the total contribution to the reaction rate from n > 6 processes is 

8 

a n (r) = 1.1) = 2.56, (124) 

n=6 

which is comparable in magnitude to 0:5(1] = 1.1) = 2.32 for the gg — ► ggg pro- 
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cesses. Therefore, they found that n-gluon processes contribute significantly to gluon 
equilibration. 

Two remarks can be made here regarding n-gluon processes. As I have dis- 
cussed earlier, the Landau-Pomeranchuk-Midgal effect is very important to regulate 
the infrared behavior of parton bremsstrahlung and to suppress soft parton radiation 
induced by multiple scatterings. It thus reduces the reaction rate of gg — > ggg. One 
would expect that this effect should also be important in multiple gluon production 
processes. Unfortunately, the analysis of the LPM effect has not been done so far for 
multiple particle radiation. In order to do so, one has to understand the space-time 
structure of the n-gluon processes. 

Secondly, the Parke- Taylor amplitude is only the tree-level amplitude for n- 
gluon processes. A complete calculation should also include all loop corrections to 
the same order in a s . Take the soft radiation limit for an illustration. In this case 
the Parke- Taylor amplitude can be factorized. Assume the nth soft gluon is radiated 
with a small angle (corresponding to a transverse momentum p?) along an energetic 
one and a fraction (z) of its momentum. The cross section for the n-gluon process 
can then be expressed in terms of that of (n — l)-gluon, 



With the gluon splitting function, P g _> gg (z) « 2N c /z for z « 1, this is much like 
the formula for successive initial and final state radiation we have encountered earlier 
in the Monte Carlo simulation [cf. Eq. (|77|) 1 except that the Sudakov form factors 
are missing, which account for the virtual corrections. One can thus include virtual 
corrections in this soft limit via multiplying the cross section by Sudakov form factors 
or introducing the '+ function' and 5-function in the splitting function. The virtual 
corrections for the complete amplitude are very complicated and so far no one has 
attempted to calculate them. However, once included, they will make the total cross 
section unitary and likely reduce the gluon multiplication rate calculated from tree 
level amplitudes. 

5 Hard Probes of the Equilibrating Parton Plasma 

So far I have discussed parton production and quark-gluon plasma formation in the 
framework of a pQCD-based model. In this framework, hard or semihard scatterings 
among partons play an essential role in liberating partons from the individual con- 
fining nucleons and driving the initially produced parton system toward equilibrium. 
The same kind of hard processes can also be used as direct probes of the early parton 




(125) 
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dynamics and the evolution of the quark-gluon plasma. Among these hard probes, 
electromagnetic signals, like direct photons and dileptons, are considered more direct 
since they can escape the dense matter without further interaction. They can thus 
reveal the dynamics of initial parton production and evolution. Similarly, open charm 
production, jet quenching due to energy loss and J/ip suppression can all provide us 
information about parton scattering and thermalization inside a parton plasma. In 
particular, preequilibrium J ftp suppression can reveal evidence of the deconfmement 
of the parton gas since the J/ip dissociation cross section with a deconfined parton 
gas is very different from that with a hadronic gas. For a recent review on signals of 
quark-gluon plasma, see Ref. ||142|| by Miiller. I will here only discuss hard probes in 
the framework of an equilibrating parton plasma. 



5.1 Open charm production 

Unlike strange quarks, charm quarks cannot be easily produced during the mixed 
and hadronic phases of the dense matter since the charm mass is much larger than 
the corresponding temperature scale. The only period when charm quarks can be 
easily produced is during the early stage of the parton evolution when the effective 
temperature is still high. At this stage, the parton gas is still not fully equilibrated 
yet so that the temperature is only an effective parameter describing the average 
momentum scale. By measuring this pre-equilibrium charm production, one can thus 
probe the initial parton density in phase space and shed light on the equilibration 
time P3 |. 



There are three different contributions to charm production in the history of 
the evolution of the parton system: (1) initial production similar to minijets; (2) 
pre-thermal production from secondary parton scatterings during the thermalization, 
r < r iso', (3) and thermal production during the parton equilibration, r > r iso . I will 
discuss them in the reversed order, starting with the charm production during the 
final stage of parton equilibration. 



5.1.1 Thermal production during equilibration 

Similar to light quark equilibration, charm quarks are produced through gluon fusion 
gg —* cc and quark antiquark annihilation qq — > cc during the evolution of the parton 
plasma. However, since the number of charm quarks is very small as compared to 
gluons and light quarks, we can neglect the back reactions, cc — > gg, qq and their 
effect on the parton evolution. Given the phase-space density of the equilibrating 
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partons, fi(k), the differential production rate is then ||144 
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(126) 



where — 2(iV 2 — 1), ^ = 2iV c , are the degeneracy factors for gluons and quarks 
(antiquarks) respectively, \A4 gg ^ c5 \ 2 , \M. qq - ¥C c\ 2 are the averaged matrix elements for 
gg — ► cc and gg — > cc processes, respectively, 
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Due to the small charm density, we can neglect the Pauli blocking of the final charm 
quarks. For large charm quark mass, M, we can approximate the phase-space density 
fi(k) by a Boltzmann distribution. We further assume that the distributions are 
boost invariant, i.e., 

fi(k) « e -*r«*%-'7) > (129) 

where 77 = 0.51n(£ + z)/(t — z) is the spatial rapidity of a space-time cells at (t,z). 
Neglecting the transverse expansion, the above assumption implies that the space- 
time cell at (£, z) have a flow velocity, u = (cosh 77, sinhr/). We can now complete the 
integral over rj in / d 4 x = ttR\ J drjdr and obtain, 
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(130) 



where Kq is the modified Bassel function and r c is the time when the temperature, 
T, drops below 200 MeV. The kinematic variables are chosen such that, 



p 2 = (M T2 cosh?/2,PT2 cos0 2 ,Pt2 sin0 2 , M T 2smhy 2 ), M T2 = \J M 2 + p\ 2 ; 

ki = k Ti (coshy ki ,cos(f) ki ,sm(t> ki ,smhy ki ), i = 1,2 . (131) 
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The center-of-mass momentum,^ = (Qt coshyg, q T , Qt sinh yg), is defined as Q = 
p + p2 = kl + k2, and 



Q 2 = s = 2 [M 2 + M T M T2 cosh(?/ - y 2 ) - PtPt2 cos <f> 2 ] , 
q% = Pt+ Pt2 + %PtPt2 cos 2 , 

Q 2 T = Q 2 + q 2 T = M| + M| 2 + 2M T M T2 cos% - y 2 ). (132) 
Using these variables and the energy-momentum conservation, we have, 
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fen = 



M T cosh(y - y fel ) + M T2 cosh(y 2 - j/ fcl ) - q T cos i? ' 
cos^i, = [p T cos (p kl + p T2 cos((p 2 - (p kl )]/q T . (133) 

In the integral over r, we shall use the time evolution of the temperature, T(r), and 
fugacities, Aj(r), as given in the previous sections. 

5.1.2 Pre-thermal production 

Before the parton distributions reach local isotropy in momentum space so that the 
rate equations can be applied to describe the equilibration of the parton system, scat- 
terings among free-streaming partons can also lead to charm production. Since the 
system during this period consists dominantly of gluons, we shall only consider gluon 
fusions. To model the phase-space distribution, we take into account the distribution 
of the initial production points which spread over a region of width 

A fe » 2 (134) 
Kt cosh y 

in z coordinate. Assuming free-streaming until r iso and neglecting the expansion in 
the transverse direction, the correlated phase-space distribution function is given by 

f(k,x) = — l — Y g(k T ,y)D(x - ztt&nhy)0(T iso coshy - t), (135) 
9GKtL A 

where 

-z 2 /2A 2 

D{x) = e -j=^-6{R A -r) (136) 
is the initial spatial distribution at t — 0, and 



71 



y/i (TeV) Y h(kr) (GeV"2) 



0.2 2.5 1754.4e- fcT /°- 9 /(A; T + 0.3) 

5.5 5.0 2.66 x 10 7 /(k T + 2.9) 6 - 4 

Table 2: Parametrizations of the momentum spectra of the initially produced partons 
in HIJING calculation. The unit of the transverse momentum k? is GeV 



is the parametrization of the parton spectrum given by HIJING simulations. The 
phase-space distribution is normalized such that lim^oo / d 3 xf(k, x)/(2n) 3 = d 3 N g /d 3 k. 
The function h(kr) and the rapidity width Y are given in Table ^|for central Au + Au 
collisions at RHIC and LHC energies. 

Substituting the phase-space distribution, into Eq. (|126[ ), and integrating over 



space and time, we obtain the charm production distribution in the pre-thermal 
period, 



dyd 2 p T 16(2n) 8 nR\ 



1 i-T-7 |2 1 

\M u 



PT2dp T 2d4>2dy 2 d(f)kidy kl —l 1 g(k T1} y k i)g(k T2 , Via) 
Jo 



2' yy ^' V2ttA._. 

t f = T iso min(coshy k i, cosh y k2 ), A tot = \J + A| 2 , (139) 
where the kinematic variables are similarly defined as in Eq. (|130 ), and in addition, 



k 2 T2 = Q? r + k^ 1 -2k T i[M T cosh(y-y kl )+M T2 cosh(y 2 -y kl )), 
cosh?/fc 2 = [Mt cosh y + Mt 2 cosh y 2 — kn cosh y k i]/kT 2 - (140) 

Note that the correlation between momentum and space-time in the phase- 
space distribution is very important as discussed by Lin and Gyulassy in Ref. ||145 



where formation time effects are also included. As we will see, this correlation will 



reduce the pre-thermal charm production as compared to the result in Ref. [143 . 



5.1.3 Initial fusion 

During the initial interaction period, charm quarks are produced together with mini- 
jets through gluon fusion and quark anti-quark annihilation. Like gluon and light 
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Figure 28: The pt distributions of initial (solid), prethermal (dot-dashed), and ther- 
mal (dashed) charm production for central Au + Au collisions at RHIC energy, 
a/s = 200 AGeV with initial conditions given in Tables p] and |2|. The dotted line 
is the thermal production assuming an initial fully equilibrated QGP at the same 
temperature. 



quark production, charm production through the initial fusion is very sensitive to the 
parton distributions inside nuclei. In addition, the cross section is also very sensitive 
to the value of the charm quark mass, M. If higher order corrections are taken into ac- 
count, the production cross section depends also on the choices of the renormalization 
and factorization scales. Detailed studies of the next-leading-order calculation ||146|| 
shows, however, that higher order corrections to the total charm production cross 
section can be accounted for by a constant K factor of about 2. This is what we will 
use next. For consistency we use M = 1.5 GeV for all calculations. Shown in Figs. ^ 
and [29| as solid lines are the initial charm production given by HIJING calculations at 
RHIC and LHC energies, with MRSD— ' |)7| parton distributions. The corresponding 
total integrated cross sections are, a c5 = 0.16 (5.75) mb at RHIC (LHC) energy, where 
nuclear shadowing of the gluon distribution function is also taken into account. In the 
HIJING calculations, higher order corrections are included via parton cascade in both 
initial and final state radiations. The resultant distributions in cc-pair momentum 
are very close to the explicit higher order calculations [|146|| . 

Plotted in Figs. |28| and |29] as dot-dashed and dashed lines are the pre-thermal 
and thermal production. In the calculation, a factor of 2 is also multiplied to the 
lowest order matrix elements of charm production. Both contributions are much 
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Figure 29: The same as in Fig. |28L except at LHC energy, y/s = 5.5 ATeV. 



smaller than the initial charm production at both energies. The pre-thermal contri- 
butions shown are also much smaller than what was found in Ref. [143] - This is be- 
cause momentum and space-time correlation was not taken into account in Ref. ||143|| 
which suppresses the pre-thermal charm production. Similar results are also found 
in a study by Lin and Gyulassy ||145 |. However, a fully equilibrated parton plasma 
(X g — X q — 1) at the same initial temperature would give an enhancement of charm 
production about 4 times higher than the initial production, shown as dotted lines 
in Figs. |2S] and 215. In this case, the enhancement not only comes from higher parton 
densities, but also from the much longer life time of the parton plasma (cf. Figs. p5| 
and |26|). 

As we have already discussed, the initial conditions in Tables [l] and |2| given 
by HIJING calculations have many uncertainties. If one increases the initial parton 
number density at RHIC energy by a factor of 4 with the same initial temperature, 
charm production from both pre-thermal and thermal sources will increase about 
a factor of 16 as shown in Fig. |30], leading to a total secondary contribution com- 
parable to the initial charm production. Alternatively, as we have discussed in the 
Section on parton evolution, we can consider the initial parton densities to be 4 times 
higher than given in Table [I] at RHIC energy but with lower initial temperature, 
To = 0.4 GeV. Accordingly, the initial phase-space distribution is also modified to: 
h(k±) = 9649.2e-^/°- 65 /0± + 0.3) from the one in Table |2|, which gives 4 times 
the initial parton density but a smaller average transverse momentum, (kj_) = 0.85 
GeV. The reduced average transverse momentum corresponds to a lower initial ef- 
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Figure 30: The p± distribution of the initial (solid), pre-thermal (dot-dashed) and 
thermal (dashed) charm production for initial parton densities 4 times higher than 
the HIJING estimate given in Table [I], with the same (ordinary lines) or reduced 
initial temperature, T = 0.4 GeV (lines with stars). 



fective temperature. This system with higher initial fugacities evolves faster toward 
equilibrium but the life-time of the deconfined phase is shorter due to the reduced 
temperature as we have discussed. The corresponding open charm production is 
shown in Fig. |30] by the lines with stars. We observe that open charm production 
from both pre-thermal and thermal contribution is reduced due to the reduction in 
initial temperature and life-time of the parton plasma, even though the initial fu- 
gacities are much higher and the evolution toward equilibrium is faster. Thus, open 
charm production is much more sensitive to the change in the initial temperature 
than to the parton fugacities. 

We note from Eqs. (|130|) and ( |138| ) that the pre-thermal and thermal charm 
production depends on the thermalization time r iso and the life time of the parton 
plasma. Therefore, by measuring the charm enhancement, we can probe both the 
initial parton phase-space distribution and the thermalization and equilibration time. 
However, in the most optimistic case, the charm production through thermal and 
pre-thermal scatterings, cannot be 40 to 50 times higher than the initial direct charm 
production, as claimed in PCM calculation ||126|| . A detailed examination exhibited 
that the copious charm production in Ref. [|126|1 is due to an overestimate of the 
intrinsic charm production. Though the intrinsic charm production is important in 
the forward direction at large Xf [147], it is strongly suppressed in the mid-rapidity 
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region due to interference among pQCD amplitudes to the same order 



I2~5l 



5.2 Dilepton and Photon Production 

The above calculation of open charm distribution does not include the effects of 
secondary scatterings suffered by the produced charm quarks. Though the secondary 
scatterings do not change the total number of charm quarks, they do modify the final 
momentum distribution. Therefore, there is no straightforward connection between 
the shape of charm distribution and the evolution history of the parton plasma. 
Dileptons and photons, on the other hand, are not subject to final state scatterings 
due to their week interaction with the medium. Their momentum distribution will 
remain intact throughout the lifetime of the system, and thus can reveal the evolution 
history. However, unlike open charm production, low-mass dileptons and small pt 
photons can be easily produced during the hadronic phase together with contributions 
from many resonance decays. Recently, it has been argued that the enhancement of 



low-mass dileptons at CERN SPS energy |148| can be interpreted by the the mass- 



shift of the p mesons ||149|| . Readers can find detailed discussions about hadronic 
contributions to dilepton production in Ref . [|150|1 . At very high dilepton mass or large 
photon px, the production is again dominated by initial hard scatterings. Therefore, 
to study dilepton and photon production from the parton phase, one has to find 
a window in phase space in which the thermal partonic contribution is dominant. 
We will find that this is not necessarily possible for all initial conditions. I will not 
discuss in detail on how to find this window and other proposed techniques [|151|| to 
single out dileptons and photons from thermal parton scatterings. Instead, I will only 
concentrate on how to calculate dilepton and photon production from an equilibrating 
parton plasma. 

Similar to open charm production, dilepton and photon production from par- 
ton scatterings can also be divided into initial, pre-thermal and thermal production. 
The thermal dilepton production rate can be written as 

= <kf££ M '?.mfMM>*), (141) 

where a q (M) = 12e 2 g Aira 2 /3M 2 is the cross section of qq — > £ + £ (M) summed over 
spin and color. Insert / d 4 PS A (P — ki — kz) and use d 4 P = (1 /2)dM 2 d 2 Prdy , we can 
obtain the differential rate as, 

f Nthc Z = Y M 2 aJM) [ d 2 P T dy 1 d^ 1 -^\fJk 1 )f,(P - h), (142) 

d 4 xdM 2 dy (2vr) 6 Y ' ^2M 2JqK JJqK h V ; 
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where k T i = M 2 /2[M T cosh(y-yi) -P T cos0i] from the requirement that (P — k 



\2 



0, and Mt = J M 2 + P^. For large values of M > T, one can approximate a boost 
invariant distribution such that 

f q (h)f q (p - h) « xyMT^Hy-WT (143) 

One can then complete the integration in momentum and space-time, d 4 x = rdrdrid 2 XT. 
The total number of thermal dileptons produced during the evolution is, similar to 
the thermal charm production in Eq. ( |130|) , 



dM 2 dy 2(2vr) 4 



r rdrM 3 TKi(M/T) J2 K a ^ M )- ( 144 ) 

J ^iso q 



This number can be calculated numerically given the time- dependence of the temper- 
ature and quark fugacity. 

Following the calculation of pre-thermal charm production, we can also com- 
pute pre-thermal dilepton yield by substituting f q (k) with a correlated phase-space 
distribution function similar to the one in Eq. (|135|) for a free-streaming system. From 
Eq. (|142|) , we can write down the dilepton yield during the pre-thermal period, 



^vrA.nt Jo 



q i 



^ e -t 2 (tanhyi-tanh?/2) 2 /2A2 o ^ ^45) 



where the time integral is similarly defines as in Eq. ( 138 ), g q = 12 is the degeneracy 
of a quark (or anti-quark) with 2 flavors, and 

k\ 2 = Pt — ^Prkn cos 4>i + k\ x 
sinh?/2 = (M T smhy — k T1 smhyi)/k T2 . (146) 

The quark distribution q(k T ,y) is obtained by multiplying the gluon distribution 
g(kT,y) in Eq. ( |137| ) by 0.07 which is determined by the the quark-to-gluon ratio in 
HIJING simulation (see Fig. 

To the lowest order, the initial dilepton production for a central AA collision 

is, 

= ? < [f ^ m2)/?(x2 ' m2) + ^ M2)f ^ M% (i4r) 
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where M is the invariant mass of the dilepton, e q the fractional charge of each quark q, 
and = e ±y Mj \fs are the fractional momenta carried by the quark and antiquark. 
Here we neglect the nuclear effects due to initial multiple interactions. One can also 
calculate the next-to-leading order (NLO) corrections which generally give a i^-factor 
about 1.1-1.4, depending on M. 

The numerical results of the above dilepton yields during three different stages 
are plotted in Figs. [H| and ^2] for central Au + Au collisions at RHIC and LHC en- 
ergies. The results of initial dilepton production are taken from a NLO calculation 
in Ref. [|152fl with MRS D— ' parton distributions. Similar to open charm production, 
we can see that dilepton production during the equilibration is much smaller than 
from initial parton scatterings due to small quark number density and the correlation 
between space and momentum for the produced partons. One also notices that ther- 
mal production at low masses becomes larger than the pre-thermal one because of 
the contribution at later time when the system is approaching chemical equilibrium 
and the quark fugacity is much larger than its initial value. We can also increase 
the initial parton density by a factor of 4 as before. The corresponding thermal and 
pre-thermal dilepton yields will also increase by roughly 16 times, which is, neverthe- 
less, still smaller than the initial production. Only in the extreme limit of complete 
thermal and chemical equilibrium initial condition can the thermal dilepton yield 
overcome the shadow of the initial Drell-Yan background in the low and intermediate 
mass range. 

This result is also consistent with earlier estimate by Eskola and Lindfors |153 |. 
My estimate here is smaller than those given by Strickland [154 and Shuryak and 
Xiong [|155|1 because of the smaller quark fugacity I used here. The effective quark 
fugacity in Ref. ||155|| is much larger because it is estimated simply from the parton dis- 
tribution function in a nucleon without taking into account smaller quark production 
cross section as compared to the gluon production. Using a parton cascade model, 
Geiger and Kapusta [ |156|| predicted a total dilepton yield which is 5 times larger than 
the underlying Drell-Yan even at M = 8 GeV region. A detailed analysis is needed 
to completely understand their results. Similarly as I discussed in multiple scattering 
and open charm production, it is very likely that the effect of finite formation time 
is also at play in this case. 

Photon production during parton equilibration can be calculated in the same 
way as done in Ref. [|154|| . One problem which is unique in photon production is that 
the quark screening mass has to be included to regulate the infrared divergency ||157|| . 
Such a regularization, however, becomes technically difficult in the calculation of pre- 
thermal photon production. However, the final result would be qualitatively similar 
to dilepton production. 
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Figure 31: The distributions of initial (solid), pre-thermal (dot-dashed), and thermal 
(dashed) dilepton production as functions of the invariant mass M for central Au+Au 
collisions at RHIC energy, y/s = 200 AGeV with initial conditions given in Tables [I] 
and 0. 
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Figure 32: The same as in Fig. |3l], except at LHC energy, ^/s = 5.5 ATeV. 
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One last note on photon and dilepton production is that background from 
hadronic decays has to be taken into account in experimental measurements. In par- 
ticular, dilepton background from charm decays is very important even at relatively 



large invariant masses ||158|| . One can even use this background, properly identified 



experimentally, to measure open charm production in pA collisions and the effect of 



nuclear shadowing [|159 



5.3 J/i/j Suppression 

I have demonstrated so far that open charm, dilepton and photon production during 
the equilibration are very sensitive to the initial parton density and the chemical com- 
position of the initially produced parton system. However, as a careful reader might 
notice, I have not addressed the question how sensitive these probes are to the color 
deconfinement of the parton plasma. In principle, charm and dilepton production in 
a hadronic gas are suppressed by form-factors as compared to a deconfined parton 
plasma at the same temperature. This is why dilepton production in the hadronic 
phase contributes mostly to the low mass region while charm production from the 
hadronic phase can be practically neglected. As we have learned, dileptons and direct 
photons are only useful when the initial parton density is high enough to produce a 
signal much larger than the background from the initial parton scatterings. Let us 
consider another process, J / ip suppression, which might be more sensitive to the color 
deconfinement but less dependent on the initial parton density. 

Using J/ip suppression as a probe of the color deconfinement |160|| requires 



that the interactions of J/ ip with hadrons and deconfined partons are different ||161| . 
Because of its small size, a heavy quarkonium can probe the short- distance properties 
of light hadrons. It is thus possible to make a parton-based calculation of the J/ip- 
hadron cross section via an operator product expansion method similar to that used in 
deeply- inelastic lepton-hadron scatterings ||161|, |162| , |163|| . The resulting J/^-hadron 



cross section can be related to the distribution function of gluons inside a hadron. The 
energy dependence of the cross section near the threshold for the break-up of a J/ip is 
determined by the large x behavior of the gluon distribution function, giving rise to 
a very small break-up cross section at low energies. Only at very high energies, this 
cross section will reach its asymptotic value of a few mb. Therefore, the dissociation 
can only occur if the gluon from the light hadron wave function is hard enough in 
the J/^'s rest frame to overcome the binding energy threshold. A hadron gas with 
temperature below 0.5 GeV cannot provide such energetic gluons to break up the 
J/ip. Therefore, a slow J/ip is very unlikely to be absorbed inside a hadron gas of 
reasonable temperature [ |161|| . 

On the other hand, a deconfined partonic system contains much harder gluons 
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Figure 33: Gluon- J/^ dissociation cross section as a function of the gluon energy q° 
in the rest frame of the J/ip. 



which can easily break up a J/ip | |161| , 164 1. The operator product expansion allows 
one to express the hadron- J /ifj inelastic cross section in terms of the convolution of 
the gluon- J/ip dissociation cross section with the gluon distribution inside the hadron 
[ 161|| . The gluon- J/ip dissociation cross section is given by [|164 



(<7%o ~ If 12 



2tt /32V 1 
Til/ m Q (e m Q )V2 (g o /eo)5 



(148) 



with q the gluon energy in the rest frame of the J/ip, vtiq the heavy quark mass, 
e = 2M D — Mj/^p the binding energy, where Mj/^ and M D are the J/tj) and D meson 



masses, respectively. As shown in Fig. [33|, the inelastic cross section shows a strong 
peak just above the break-up threshold of the gluon energy, with a maximum value 
of about 3 mb. Low-momentum gluons have neither the resolution to distinguish the 
heavy constituent quarks nor the energy to excite them to the continuum. 

In the pre-equilibrium stage, i.e., before the partons have reached equilibrium, 
the average parton transverse momentum is sufficiently large to break up a J/ip, pro- 
vided the partons are deconfined. The dissociation of J/ip will continue during the 
whole equilibration process until the effective temperature drops below a certain value 
or the beginning of hadronization, whichever takes place first. Therefore measure- 
ments of J/ip suppression can probe the deconfinement of the early partonic system 
and shed light on the subsequent equilibration process, provided that possible nuclear 
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effects on the production of QQ pairs and on pre-resonance charmonium states are 
understood and taken into account. 

Let us consider J /if> suppression in the central rapidity region (yj/$ — 0). In 
this case, the J/ if) will move in the transverse direction with a four- velocity 



u 



(M T ,P T ,0)/Mj/^ (149) 



where Mt = \J Pt + ^j/»/> * s defined as the J/^'s transverse mass. A gluon with a 
four-momentum k = (k°, k) in the rest frame of the parton gas has an energy q° = k-u 
in the rest frame of the J /if). Using Eq. ( |148[ ), one can find out the thermal averaged 
gluon- J /if) dissociation cross section ID.65I1 , 

3 C(3) PtM t T 6 m Q 

00 r oo ( r _lW2 _ 

VtJ dx^ J_( e -a n x _ e -a+xj ( 15Q ) 

n=l h X 

with T n = T/n and 

at = (151) 

J-n JMj/^ 

Using the thermal cross section just obtained, we can now calculate the survival 
probability of J /if) in an equilibrating parton plasma, Again, we will neglect the 
transverse expansion and consider only longitudinal expansion. A J /if) produced at 
point r with velocity v in the transverse direction will travel a distance 



rf = _rcos0+ v / J R2 i _ r 2( 1 _ cos 2 0) ( 152 ) 

in the time interval = MtcI/Pt before it escapes from a gluon gas of transverse 
extension Ra] here cos0 = v ■ r. Suppose the system evolves in a deconfined state 
until the temperature drops below a certain value, which we assume to be 200 MeV. 
The total amount of time the J /if) remains inside a deconfined parton gas is the 
smaller one of the two times and tf, the life-time of the parton gas. Assume 
that the initial production rate of the J /if) is proportional to the number of binary 
nucleon-nucleon interactions at impact-parameter r, N^(r) = A 2 (l — r 2 / R\)/2ttR\. 
The survival probability of the J / if) averaged over its initial position and direction in 
an equilibrating parton gas is 

q(p , I d 2 r(R 2 A ~ r 2 ) exp [- dTn g (r)(v Tel a(k ■ u)) k ] 

S{Pt) ~ Jd 2 r(R A -r 2 ) ' (153) 
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Figure 34: (a) The survival probability of J/ip in an equilibrating parton plasma at 
RHIC and LHC energies with initial conditions given by HIJING simulation in [I] (b) 
and for an initially equilibrated plasma at the same temperatures. 

where 

t min = min(^,t/) , (154) 

and n g {r) is the gluon number density at a given time r. With the time dependence 
of the gluon fugacity and temperature we can then evaluate numerically the survival 
probability of a J/ip in an equilibrating parton plasma. 

Shown in Fig. |34"| a are the J/ip survival probabilities in the deconfmed and 
equilibrating parton plasma at RHIC and LHC energies. We find that there is stronger 
J / ip suppression at LHC than at RHIC energy, due both to the higher initial parton 
densities and longer life-time of the parton plasma. The increase of the survival 
probabilities with J/ip's transverse momentum is a consequence of the decrease of 
the thermal cross section with increasing Pt at high temperatures, and the shorter 
time spent by a higher- Pt J / ip inside the parton plasma, an effect first considered in 
Ref.jT66|. 



As in open charm and dilepton production, we multiply initial parton number 
densities by a factor of 4, thus increasing initial parton fugacities to test the sensitivity 
of J/ip suppression to the uncertainties in the initial conditions. The corresponding 



survival probabilities, shown in Fig. EB as solid lines, are much lower. If the un- 



certainties in initial conditions are caused by the soft parton production from the 
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Figure 35: The survival probability of J /if) in an equilibrating parton plasma with 
initial parton densities 4 times higher than HIJING estimate given in |I| but with the 
same (solid line), or reduced initial temperature, T = 0.4 GeV (dashed line). 



color mean fields, the initial effective temperature will decrease. Therefore, we can 
alternatively increase the initial parton density by a factor of 4 and at the same time 
decrease T to 0.4 and 0.72 GeV at RHIC and LHC energies, respectively. This leads 
also to higher initial fugacities. Comparing the solid and dashed lines shows that 
the J /if) suppression is less sensitive to the variation of the initial temperature and 
fugacities as long as the parton densities are fixed. In the most extreme case of a fully 
equilibrated parton plasma (unit initial fugacities), J /if) is much more suppressed, as 
shown in Fig. |3"^}b, because of both the higher parton density and the longer life time 
of the system. 

At high temperature, the thermal-averaged gluon-J/^ dissociation cross sec- 
tion decreases monotonically with Pt since the averaged collision energy is already 
above the break-up threshold. One can translate this into a J /if) survival probability 
monotonically increasing with Pt. However, at low temperature blow 300 MeV, the 
Pt dependence of the thermally averaged gluon-J/^ dissociation cross section also 
displays the threshold peak as in Fig. |33]. The peak moves to larger values of Pt 
at a lower temperature ||165|| . Therefore, due to the large values of the cross section 
around the threshold peak, the resultant survival probability integrated over the evo- 
lution time at low temperatures is quite flat at small Pt- The final shape of the Pt 
dependence of the survival probability at small Pt depends on the entire evolution 
history of the parton system, especially the relative lengths of time the system spends 
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in the high and low temperature stages. The smaller Pt dependence at small Pt seen 
in Figs. [£| and ^ is therefore due to the J /if) dissociation during the late stage of 
the evolution when the temperature is below 300 MeV [ 165|| . For a parton system 
with a low initial temperature (below 300 MeV), the Pt dependence of the survival 
probability should be even flatter. One can therefore use the Pt dependence to shed 
light on the initial temperature and the evolution history of the system. 

We have seen that the large average momentum in the hot gluon gas enables 
gluons to break up the J /if), while hadron matter at reasonable temperature does 
not provide sufficiently hard gluons. A substantial J /if) suppression occurs in such 
a non-equilibrium partonic medium, though smaller than that in a fully equilibrated 
parton plasma. I should emphasize that, in addition to the J/if> dissociation during 
the equilibration of the parton plasma, there are other possible sources of suppression 
for the actually observed J/ip's, such as nuclear modifications of the QQ production 
process e.g., through modified gluon distributions in a nucleus p5| , |100| , multiple scat- 



tering accompanied by energy loss [ 167 , or a suppression of the nascent J /if) before it 
forms an actual physical resonance ||168| , |169|| . Such effects would cause J /if) suppres- 
sion in addition to what we have obtained from the equilibrating parton plasma and 
modify the transverse momentum dependence of J/if> suppression ||170| . Moreover, 
there should also be J /if) pre-thermal suppression. On the other hand, gluon fusion 
could also result in J / if) production during the evolution of the parton system, similar 
to the enhancement of open charm. A consistent study of J /if) suppression should 
include the pre-equilibrium production in a form of a master rate equation. 



5.4 Jet Quenching and Monojet Production 

Complementary to J/if> suppression, a study of high pt jets and their propagation 
inside the medium can also probe the structure of the dense matter and possibly the 
phase transition, since these high pt jets are produced on a very short time scale 
as compared to the soft processes and their production rates and spectrum can be 
reliably calculated via pQCD. For example, an enhanced acoplanarity of two back-to- 
back jets in nuclear collisions can be used to study multiple scatterings of a parton 
inside a dense medium 
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What jets could also probe in high energy heavy ion 
collisions is the stopping power, dE/dz, of the dense matter for high energy quarks 
and gluons |88|, |108|, |109||. That stopping power in turn is controlled by the color 



screening mass fiD in that medium. A possible rapid change of \xb near the phase 
transition point could lead to a variation of jet quenching phenomena which may serve 



as signatures of QGP formation [ 109|| . The energy loss, dE/dz, of partons through 
interaction is also closely related to the thermalization and equilibration of partonic 
system as we have discussed. 
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5.4.1 Energy loss of a fast parton in QGP 

When an energetic parton propagates through a quark-gluon plasma, it suffers both 
elastic energy loss, dE e \/dz due to simple scatterings, and radiative energy loss 
dE ra d/dz due to induced radiations from multiple scatterings. Studies [|172| , |173|1 



of the elastic energy loss at finite temperatures give a logarithmic energy depen- 
dence of dE e \/dz. Its value is also quite small, for example, 0.2 — 0.5 GeV/fm for 
a quark jet with E = 30 GeV. The dominant source of energy loss is expected to 
be caused by induced radiations. The estimates of dE ra _ d /dz in the past vary widely 
174], [167| , |175| , |176|| from energy independent to dE/dz oc E. M. Gyulassy and I 



have analyzed soft gluon bremsstrahlung |27], [3T| and found that interference due to 
multiple scatterings is very important to calculate the radiative energy loss. Similarly 
to the parton equilibration rate, the LPM effect, i.e., the suppression of soft gluons 
whose formation times are much larger than the mean free path of parton scatterings, 
also reduces the radiative energy loss inside a dense QCD medium. Dokshitzer et al. 
[ ft2| have made a more complete analysis of the QCD radiation amplitude, including 
also the rescatterings of the radiated gluon. Their results also shown suppression 
of soft gluon radiation but with a square-root energy-dependence of the energy loss. 
Their energy loss, dE/dz, is also proportional to the distance z the propagating par- 
ton has traveled. To demonstrate the LPM effect on energy loss and jet quenching, 
we will use a simple estimate in this paper. 

To estimate the radiative energy loss, we can use the regularized gluon spec- 
trum induced by a single scattering in Eq. (|1 1 1|) . We then again take into account 
the LPM effect in QCD due to color interference by approximating the suppression 
factor with a step function, 0(\f — tqcd)- This step function restricts the phase space 
of the radiation to a region in which the effective formation time tqcd(^) must be 
smaller than the mean free path A/. Radiations not fulfilling this requirement will be 
suppressed by the destructive interference. The radiative energy loss suffered by the 
jet parton per collision is then, 

f dn 

A£ rad = J d 2 k ± dy d2k ^ dy k± coshy9(Xf - r QCD (k))9(E - k ± coshy), (155) 

where tqcd(^) is given by Eq. ( |110| ), E is the jet parton energy, and the regular- 
ized gluon density distribution, dn g /dyd 2 k±, induced by a single scattering is given 
by Eq. ( |111| ). Since the transverse momentum transfer qj_ is the result of elastic 
scatterings, we have to average any function /(q_i_) of q_i_ by the elastic cross section, 

(f^±)) = ~ r%l^/(q±), (156) 
Oi Jii 2 D aq ± 



N6 



where s ~ 6ET is the average cm. energy squared for the scattering of a jet parton 
with energy E off the thermal partons at temperature T. Since we now only consider 
a fully equilibrated quark-gluon plasma, the Debye screening mass will be given by 
Eq. (|108|) with X g = 1, i.e., \j? d = Ana s T 2 . For the dominant small angle scattering, 
the elastic cross sections are, 

(157) 

where C % el = 9/4, 1, 4/9 respectively for gg, gq and qq scatterings. This average can 
be approximated by replacing qj_ in the numerator of Eq. ( |111[ ) with its average value, 



(g±> = (158) 

In the denominator, we simply replace by y? D after the angular integration. The 
remaining integration in Eq. (|155|) over the restricted phase space approximately leads 
to the simple analytic formula, 



AE rad « ^(gi) ^Jx + -^J 2 , (159) 
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(161) 



In small k± regime, the phase space is mainly restricted by a small effective forma- 
tion time, tqcd < A/, which gives the first term proportional to A/. For large k±, the 
radiation becomes additive in a restricted phase space constrained by energy conser- 
vation. That region contributes to the second term which appears to be proportional 
to the incident energy E. However, in the high energy limit the function X 2 oc l/E 
and hence the radiated energy loss grows only as log 2 E. 

The above derivation assumed that the mean free path is much larger than 
the interaction range specified by 1//xd- As we shall discuss below, this is satisfied in 
a quark gluon plasma at least in the weak coupling limit. Therefore, we can neglect 
the second term in T x . For a high energy jet parton, E 3> fin, we can also neglect 
the second term in X 2 . The resulting radiative energy loss reduces in that case to the 
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simple form, 



dE, 



rad 



AE- 



rad 



Co a.. 



dz 



7T 




(162) 



which depends on a dimensionless variable, 



(163) 



We see that the radiative energy loss dE rad /dz thus obtained interpolates be- 
tween the factorization and Bethe-Heitler limits as a function of the dimensionless 
variable £. In the factorization limit, we fix //dA/ ^> 1 and let — > oo, so that £ 3> 1. 
In this case, we can neglect the second term in Eq. (|162j) and have, 



dR 



rad 



dz 



C 2 a s 2 ( 2r 2 E \ 



7T 



£»1. 



(164) 



Thus, the radiative energy loss in the factorization limit has only a logarithmic energy 
dependence (in addition to the energy dependence of (q±))- Due to the non-Abelian 
nature of the color interference, the resultant energy loss for a gluon (C 2 = Ca) is 
9/4 times larger than that for a quark (C 2 = Cp). In the other extreme limit, we 
fix E and let //dA/ — ► oo, so that £ <C 1. In this case, the mean-free-path exceeds 
the effective formation time. The radiation from each scattering adds up. We then 
recover the linear dependence of the energy loss dE Tad /dz on the incident energy E 
(modulo logarithms), 
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dz 27rA/ \j? d 
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'2^ D X f 
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(165) 



as in the Bethe-Heitler formula. In both cases, the radiative energy loss is proportional 
to the average of the transverse momentum transfer, (qj_), which is controlled by the 
color screening mass as in Eq. ( |158| ). 

To see more clearly how the factorization limit is approached, we now estimate 
£ for a parton propagating inside a high temperature quark-gluon plasma. From 
Eq. (|157|) and the quark and gluon densities in an ideal system at a temperature T, 
the mean free path (for 3 quark flavors) is 



1/ = a qqPq + a qgPg 



2-kq 2 T 3 
^4 x 7C(3)^, 



7T Z 



(166) 
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Figure 36: The energy dependence of energy loss, dE/dz, of a quark with energy E 
inside a quark-gluon plasma at temperature T = 300 MeV. A weak coupling a s = 0.3 
is used. The solid line is the full expression and the dashed line is the factorization 
limit of the radiative energy loss. The dot-dashed line is the elastic energy loss. 



, , 27ra 2 T 3 

W = v q9 p q + cr 9gPg « -^9 x 7C(3) — , (167) 

where C(3) w 1.2. 

Using Eqs. (|166|) and ( |167|) and the perturbative color electric screening mass, 



fijj = 4ira s T 2 , we see that £ appearing in the logarithms has a common energy and 
temperature dependence for both quarks and gluons, 

£ = , 2 o = -—^-t)-— ~ —7 — , for both q and o. (168) 
S X f li 2 D 16vr 3 T 2vr 3 T' y y v ; 

With the above expression for £, we plot the radiative energy loss in Fig. ^ 
as a function of the jet energy inside a plasma at temperature T = 300 MeV with 
a s = 0.3. The solid line is the full expression in Eq. ( |162| ) while the dashed line is the 
factorization limit corresponding to the first term in Eq. ( |162| ). We see that Eq. ( |164| ) 



approximates Eq. ( |162j ) quite well in this parameter range. The energy dependence 
of the radiative energy loss is due to the double logarithmic function in the formula 
one of which comes from the energy dependence of the average transverse momentum 
(ql) in Eq. (|T58|). 

The energy loss of a quark in dense matter due to elastic scattering was first 



estimated by Bjorken ||172]| and later was studied in detail ||1 73|| in terms of finite 
temperature QCD. For our purpose, a simple estimate taking into account both the 
thermal average and color screening will suffice. In terms of elastic cross sections and 
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the density distributions for quarks and gluons in a plasma, we have, 



"4 d ^dl Pi " = {q±)ai{ ^ h (169) 

where v « q\/2u is the energy transfer of the jet parton to a thermal parton with en- 
ergy uj during an elastic scattering, (q\) is the average transverse momentum transfer 
given by Eq. ( |158| ). Similar to Eqs. ( |166| ) and ( |167| ), we have, 



•qq 



The elastic energy loss of a fast parton inside a quark gluon plasma at a temperature 
T is then given by, 



For comparison, we plot this elastic energy loss in Fig. |36|. In general, it is much 
smaller than the radiative energy loss and has a weaker energy dependence (single 
logarithmic). 



Using Eqs. ( |158| ), ( |164j) and ( |168j ), the total energy loss can be expressed as, 



dE = dE^ + dE^ _ 2 ^ 3ET / 9E_ 3n 2 a 

dz dz dz 7i ^ D 2p 2 D I 7r 3 T 2p? D 



It is interesting to note that both the elastic and radiative energy loss have the same 
color coefficient C^- For high energy partons, the radiative energy loss dominates over 
the elastic one. For E = 30 GeV, T = 300 MeV, and a s ~ 0.3, the total energy loss 
for a propagating quark is dE/dz w 3.6 GeV/fm. Only about 25% of this amount 
comes from elastic energy loss. 

Comparing to the phenomenological string tension k ~ 1 — 2 GeV/fm, which 
can be considered as the dE/dz of quarks in the nonperturbative QCD vacuum ||178 



one cannot really distinguish jet quenching caused by energy loss in a quark gluon 
plasma from that in a hadronic matter. However, since the energy loss is very sensitive 
to the infrared cutoff scale p£,(T), the study of jet quenching may be used as a 
signature of a quark gluon plasma if the Pd(T) dependence of dE/dz can still be 
given by Eq. ( |173j ) and if there is a significant variation in po{T) near the phase 
transition temperature T c . Up to now, little is known about Pd{T) for physical QCD. 
However, lattice simulations [ p.77[ | of a pure gauge QCD indicate a rapid decrease 
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Figure 37: Lattice result of the Debye screen mass JL77J, //p, as a function of tem- 
perature. The curve is a parametrization used in the following calculation. 



of Hd(T) near T c as shown in Fig. |37], where the weak coupling limit, /x(T) = ^T, 
is indeed recovered at high temperatures. Since a sharp reduction of Hd(T) would 
effectively suppress the energy loss according to our formula in Eq. ( |173| ), a long 
duration of the mixed phase in a near-first order phase transition could lead to a 
decrease of jet quenching. 

In the following, we will use a parametrization of the lattice result as a phe- 
nomenological guide to describe the temperature dependence of the color screening 
mass Hd(T) both near and above T c . This will give us the temperature dependence of 
the energy loss according to Eq. (|173|) . The coupling constant which corresponds to 
/Ad(T) = gT in Fig. [37| at high temperatures is a s ~ 0.3. We shall freeze the coupling 
constant to this value at all temperatures in calculating dE/dz in the QGP phase. 
In hadronic phase at T = 0, we assume a constant energy loss given by the string 
tension, k ~ 1 GeV/fm, which is consistent with the deep inelastic lepton-nucleus 
scattering data ||179|| . To extrapolate from T = to T c , one may choose a smooth 
form, dE/dz 
jet with E = 



-> 1 GeV/fm as T — ► 0. This scenario is plotted in Fig. 38 for a quark 
30 GeV as set 1, in which there is a dip in dE/dz during the phase 
transition due to the small value of \lb at T c . We have assumed a QCD phase tran- 
sition temperature T c = 0.188 GeV. Alternatively, one may consider a scenario (set 
2 in Fig. [38|) in which dE/dz = cont. = l GeV/fm in the hadronic phase, i.e., there 
is a discontinuity in dE/dz at T = T c . For comparison, we also study the case that 
dE/dz in the T = limit is 2 GeV/fm (set 3). 
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Figure 38: Three scenarios of energy loss as functions of temperature. 



5.4.2 Suppression of dijets 

We now discuss the implications of the temperature dependence of dE/dz on the 
suppression of di-jets and on the production rate of mono-jets in ultrarelativistic 
nuclear collisions. Similarly to Eq. (|19D, the initial rate of back-to-back parton jets 
of transverse momentum px and rapidities y\ and y 2 produced in hard scattering 
processes in central AA collisions is, 

dijet ' ' r >- ~^ '- » " * — ■* ~ » *> > » ' o \ tiu 



., f d 2 bT\(b)K xifa/N(xi,p T )x 2 fb/N(x2,P T ) -^-j? — (s,i,u) (174) 

dp T dy x dy 2 J a ^ cd dt 

where we have neglected the nuclear shadowing effect for large pr jet production, i.e., 
f a /A(x,Q 2 ) ~ fa/N(%,Q 2 )TA(b), T^ip) is the nuclear thickness function as defined 
in Eq. (|66D. The integration is performed over the plane transverse to the beam 
direction, b is the transverse coordinate of the production point in the transverse 
plane. 

In the presence of dense matter, the high energy partons produced by a hard 
scattering process can interact in the medium and thereby lose energy which will not 
show up in the final state of the jet but contribute to the soft background. Thus, 
the number of QCD hard jets to be observed in AA collisions will be reduced as 
compared to the expression Eq. (|174 ) (in Ref. ||108| the term "jet quenching" has 



been introduced for this effect). Taking into account interactions and energy loss in 
the medium, the rate of di-jets of rapidities yi,y 2 and transverse momenta Pt\,Vti 
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reads 



dp T1 dp T 2dyidy 2 



d 2 b J dp 2 T Tl(b)K ^ x 1 f a/N (x 1 ,p 2 r ) x 2 f b/N {x 2 ,p 2 T ) 

a,b,c,d 

(s, i, u)S(pri ~Pt + AE Ct i(b, p T , yi)) 



dt 

■ S(p T 2 ~Pt + AE dt2 (b, p T , y 2 )) (175) 

The energy losses AE C ^ and AE^ depend not only on the momenta and on the 
production point of the partons c and d but also on the space-time evolution of the 
dense matter. If we limit ourselves to jets in the central region, y\ = y 2 = 0, the total 
energy loss of a massless parton inside the dense matter is 

AE = Hdr^ir) (176) 
J to dz 

where r is the formation time of the dense matter and Tf is the proper time when 
the parton escapes from the dense medium, Tf = (R 2 — b 2 sin 2 if)) 1 / 2 — b cos ip, where 
■0 is the angle between b and pt- 

To evaluate the integrated energy loss suffered by a fast parton, we now need to 
know the the time dependence of the temperature. We shall assume a one-dimensional 
scaling expansion of the system ||118||, thus the time dependence of the temperature 



in both phases is given by T ~ 1/r 1 / 3 . To determine the duration of the mixed phase 
of a first order phase transition in which the temperature is constant, T c , we adopt a 
bag model equation of state (EOS) for the quark-gluon-plasma (QGP) phase: 

37 1 
e(T) = jjV T 4 + B P(T) = -[e(T) - AB] (177) 

where e, P and B are the energy density, the pressure and the bag constant, respec- 
tively. Define the energy densities at the critical temperature in the QGP phase and 
in the hadronic phase, Eq and eh, respectively. The system then enters the mixed 
phase at the proper time tq = t Tq/T c 3 , and goes from the mixed phase into the 
hadronic phase at t h = {eq + P c )tq/(eh + P c ), where P c = (eq — 4B)/3 is the pres- 
sure at T c . In the following, we shall use Eq =2.5 GeV/fm 3 , e h =0.5 GeV/fm 3 , and 
B =0.5 GeV/fm 3 , which implies a critical temperature T c =0.188 GeV. 

In order to achieve maximum sensitivity to the quenching effect it is useful to 
consider the total transverse energy E tot = p T \ + Pt2 for a pair of back-to-back jets. 
The quenching ratio is defined as the number of dijets of energy Etot as measured in 
the experiment, divided by the corresponding number expected in the absence of jet 
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Figure 39: The suppression rate as a function of the initial energy density e$ for a 
dijet with the total energy E T = 30, 60 GeV in a central U + U collisions at y/s = 200 
GeV. 
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Figure 40: The same as Fig. EM except at y/s = 6 TeV 
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interaction with the dense medium, 



Raa{Eu 



J dpri J dp T 2 5(E tot -prx - PT2) 



dN, 



dijet 



dp T1 dp T2 dy 1 dy 2 



J dp 2 T 5{E tot - 2p T ) 



dN (0) 

ttJV dijet 



dp^dy x dy 2 



2/1=2/2=0 



for dijets at y\ = y 2 = in the CM frame, 

Figs. ^ and ^ show the dependence of the quenching rate on the initial 
energy density e for central U + U collisions at yfs =200 AGeV and y/s =6 ATeV, 
respectively, for given total dijet energies. The results were obtained by evaluating 
Eqs. ( |174j ) and (|175[ ) for the Duke-Owens structure functions set 1 |55| and the EOS 
discussed above. The initial proper time was taken to be r = 1/Tq. The non- 
trivial dependence of RAA(E tot ) on Eo reflects the temperature dependence of dE/dz 
as shown in Fig. For small initial energy densities, the high energy partons of 
the dijet suffer a temperature-independent energy loss per unit distance, dE/dz = 1 
GeV/fm, in the hadronic phase of the dense matter at temperatures well below T c . 
As Eo increases the partons will spend more and more time in the mixed phase. For 
all scenarios , dE/dz takes its minimum value in the mixed phase and consequently, 
the quenching rate increases with increasing initial energy density, i.e., one observes 
less quenching. On the other hand, in the absence of such a local minimum, no 
increase of the dijet rate with e q occurs. As Eq further increases, the leading partons 
of the jet probe the QGP phase at higher temperatures. Since dE/dz grows with 
T for T > T c , RAA(Etot) decreases with e in this region. All three curves come 
together asymptotically as they correspond to the same dE/dz in the QGP phase. 
Thus a maximum of RAA{Etot) a t intermediate Eo reflects a minimum of dE/dz at 
the critical temperature of the phase transition. 



5.4.3 Monojet production 

Depending on the position of the initial hard parton scattering and the direction of 
their momenta after the interaction, it can happen that one of the high energy partons 
has to traverse a large distance in the dense matter while the other one travels only a 
short distance in the medium. If the first parton loses such a large amount of energy 
that the corresponding jet can no longer be distinguished from the soft and semi-hard 
background, one observes a mono-jet. To be specific, let us define a mono-jet as a 
dijet where one of the jets has a transverse momentum below a given value PTcut- 
The number of mono-jets of transverse momentum and rapidity y per central AA 
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Figure 41: Monojet fraction as a function of the initial energy density Eq with a fixed 
Pt = 30 GeV and different values of pTcut for central U + U collisions at y/s = 200 
GeV, and 6 TeV. 



collision at y = is given by [cf. Eq. ( |T75| )1 
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(179) 



where one has to integrate over the (experimentally unknown) rapidity and transverse 
momentum of the second parton. 

For large Eq, the number of mono-jets must decrease because of the increase 
of dE/dz with temperature in the QGP phase which leads to the suppression of all 
jets. In this context it is useful to define the mono-jet fraction, i.e., the mono-jet rate 
divided by the total jet rate, at a given transverse momentum p T and rapidity y = 0, 



Rmono{PT,PT2 < PTcuuD = 0) 



dy 2 / dp T 2 6(prcut -Pt%) 



dN Ai]et 



dp T dp T2 dydy 2 



y=0 



dy2 / dp T2 



dN d ; 



dijet 



dp T dp T2 dydy 2 



y=o 



(180) 

In Fig. [4T], -R mon o(PT = 30 GeV/c,px2 < PTcnt, y = 0) is plotted as a function of e for 
scenario 1. Clearly, the unquenching of jets associated with the minimum of dE/dz 
at the critical temperature leads to a dip in the mono-jet rate which corresponds to 
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the peak in the di-jet quenching rate (cf. Figs. [39| and EDI) . In the limit of large initial 
energy densities -Rmono — ► 1, i-e., all jets are mono-jets at large Eq. The reason is 
that eventually dE/dz becomes so large that, except for negligible fraction of parton 
scatterings whose production points are close to the nuclear surface in the transverse 
plane, only one of the two jets can survive. 

Since the jet production rates depend on energy, A which determine the trans- 
verse size of the dense matter, Ra, and the initial energy density 6q, one would like to 
fixed the first two variables and isolate only the dependence on e - m experiments, 
there is no practical method to vary the initial energy density Eq for both fixed yfs 
and A. In order to best study jet quenching and monojet production as a probe of 
the phase transition, one can vary the initial energy density by either changing s/s 
for fixed beams or changing A number for fixed energy ||110| . 



6 Summary and Discussion 

In this report, I presented a pQCD-based picture of ultrarelativistic heavy-ion col- 
lisions. In this framework, a nucleus in the infinite momentum frame consists of 
many partons (quarks and gluons). The interactions among these partons can be 
divided into perturbative, which can be described by pQCD calculation, and non- 
perturbative, which can only be modeled phenomenologically. I have demonstrated 
that pQCD processes dominate the underlying dynamics of heavy-ion collisions at 
extremely high energies. I argued that the soft component of a Pomeron exchange 
would be suppressed in the presence of a dense and hot partonic plasma, though it is 
still important in the initial parton scatterings in the early stage of the heavy ion col- 
lisions and is responsible for the initial soft parton production. It is then reasonable 
to assume that the evolution of the initially produced partons can be described by 
pQCD processes. Using the initial conditions estimated by the HIJING Monte Carlo 
model, the following picture emerges: 

(1) During the early stage in ultrarelativistic heavy- ion collisions, hard or 
semihard parton scatterings, which happen in a time scale of about 0.2 fm/c, produce 
a hot and undersaturated parton gas. This parton gas is dominated by gluons and is 
far from chemical equilibrium. Multiple hard scatterings suffered by a single parton 
during this short period of time when the beam partons pass through each other 
are suppressed due to the interference embedded in the Glauber formula for multiple 
scatterings. This leads to the predicted disapparence of the Cronin effect. Interference 
and parton fusion also lead to the depletion of small x partons in the effective parton 
distributions inside a nucleus. This nuclear shadowing of parton distributions reduces 
the initial parton production. 
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(2) After the two beams of partons pass through each other, the produced 
parton gas in the central rapidity region starts its evolution toward (kinetic) ther- 
malization and (chemical) equilibration through elastic scatterings and induced radi- 
ations. The kinematic separation of partons in the central slab of about 1 fm through 
free-streaming gives an estimate of the time scale r iso ~ 0.5 — 0.7 fm/c, when lo- 
cal isotropy in momentum distributions is reached. Further evolution of the parton 
gas toward a fully equilibrated parton plasma is dictated by the parton proliferation 
through induced radiation and gluon fusion. Though the gluon equilibration rate is 
reduced by the inclusion of the Landau-Pomeranchuk-Migdal effect, gluon fugacity 
still increases rapidly toward its equilibrium value. Due to the consumption of energy 
by the additional parton production, the effective temperature of the parton plasma 
cools down considerably faster than the ideal Bjorken's scaling solution. Therefore, 
the life time of the plasma is reduced to 4 - 6 fm/c before the temperature drops 
below the QCD phase transition temperature. 

(3) The evolution of the quark distribution always lags behind that of gluons 
due to a smaller equilibration rate and the initial density. For heavy quarks, the 
equilibration rate is even smaller. Take charm quarks for example. The thermal 
production during the equilibration period is much smaller than the initial direct 
production, due to the small initial gluon fugacity and the short life time during 
which the temperature remains high enough to produce charm quarks. For the same 
reason, dilepton and photon production during the evolution of the parton plasma 
is also small because of small quark number density. Therefore, observation of large 
charm and dilepton enhancement would imply high initial gluon and quark density 
and thus a longer life time of the parton plasma [|143j , |144j , |145| . 

(4) Even though the initial parton system is not in a full equilibrium, a study of 
color screening ||135fl shows that the system is already in a deconfmed state with large 
average momentum (or effective temperature). Such a deconfmed parton system, 
though not in equilibrium, will dissociate hadronic states like J/ip. It was shown that 
J I ip can be substantially suppressed during the evolution of the parton plasma toward 
equilibrium. The measurements of this suppression can reveal the initial conditions 
and the evolution history of the parton plasma in high-energy heavy-ion collisions. 

Readers should be reminded that the above picture is the result of only a 
qualitative estimate based on a pQCD-based model. Many uncertainties prevent us 
from making a more quantitative calculation. All in all, these uncertainties arise from 
our ignorance of the nonperturbative physics and our inability to calculate the soft 
processes in the framework of QCD. In the pQCD-based model reviewed here, the un- 
certainties really lie in the cut-off, p , which is supposed to separate nonperturbative 
soft interactions from perturbative hard processes. Since soft and hard physics do 
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not have a definite boundary, the resultant parton production from hard or semihard 
processes is very sensitive to the cut-off. The accompanying soft parton production 
is not known in this model and may only be estimated by effective models like the 
color flux-tube model ||61|| . Recently, it has been argued that the presence of larger 
partons can provide screening to the parton production with smaller p? ||180|| . The 
screening mass can then replace the ad hoc cut-off and regulates the parton produc- 
tion cross section self-consistently. There have been recent developments in the field 
theory of particle production from mean fields |24] . Such a field theoretical approach 
to particle production could be the ultimate and consistent way to address the pro- 
duction and formation of a parton plasma in heavy-ion collisions, once the collision 
terms are included and the theory is generalized to the QCD case with rather random 
color fields. 

The other nonperturbative uncertainty comes from the parton distribution 
functions inside a nucleus. Though quark distributions have been measured in deeply 
inelastic lepton-nucleon (nucleus) collisions, the gluon distribution is not known very 
well. Furthermore, nuclear modification of the gluon distribution has not at all been 
measured experimentally. Since gluon scatterings are the dominant QCD processes 
in parton production, our results are very sensitive to the gluon distribution. On 
top of that, nuclear shadowing will cause additional 20-50% uncertainty in the initial 
parton density. One way to reduce this uncertainty is to measure pp and pA collisions 
systematically as proposed in Ref. |88], |159 . 



For a thermalized parton plasma, the uncertainties related to nonperturbative 
physics will be reduced, since the soft processes (long distance interactions) will be 
screened by the interactions with the hot medium. For example, the original diver- 
gent cross sections are regularized by the color screening mass fin and the infrared 
and collinear singularities associated with gluon radiations can be regularized by the 
LPM suppression. Recent developments in QCD transport theory p5| have systemat- 
ically studied these problems and formally included the resummation of hard thermal 
loops, first proposed by Braaten and Pisarski [133] in finite temperature field theory. 
However, one finds that, even though the electric sector of the strong interaction is 
regularized by an effective Debye screening mass, the magnetic interaction is still not 
regularized. For most of the transport coefficients, this is not a problem since the 
cancellation between in and out-state scatterings makes the effective cross section less 
divergent and the imaginary part of the gluon self-energy is enough to regulate the 
infrared divergency | 132 , |130| , 122 |. This is so-called dynamic screening by Landau 
damping. For some other quantities, like gluon damping rate ||181|| and color diffusion 
rate | |131| , [132| |, one still has to invoke a magnetic screening mass which will always 
be a nonperturbative quantity. Higher order corrections to the Debye screening mass 



99 



also depend on the non-perturbative magnetic screening mass [182] . Therefore, some 
aspect of the parton thermalization and equilibration are still influenced by nonper- 
turbative physics. 

Throughout the production and evolution of the partonic system, interference 
effects play an important role in multiple collision processes, the correct implementa- 
tion of which will be a great challenge to any transport model based on classical parton 
cascades. First, the interference between different amplitudes in multiple parton scat- 
terings leads to the Glauber formula which suppresses multiple hard scatterings of 
beam partons. Second, the LPM effect, i.e., the destructive interference among dif- 
ferent amplitudes of gluon radiation induced by multiple scatterings, suppresses soft 
gluon radiation whose effective formation time, tqcd, is larger than the mean free 
path, A/, of the parton scatterings. This leads both to reduced gluon equilibration 
rates and reduced radiative energy loss. A detailed analysis of the radiation ampli- 
tude has revealed the underlying physics which contradicts the intuitive picture of the 
classical parton cascade model. A propagating parton can no longer be considered 
as time-like with a degrading virtuality as it radiates between multiple scatterings, 
since the parton can be space-like (q 2 < for initial state radiation) in one amplitude 
and time-like (q 2 > for final state radiation) in another. It is the cancellation be- 
tween these amplitudes that suppresses the soft radiation with long formation time. 
One way to incorporate the LPM effect in the parton interaction simulation, as I 
have illustrated, is to consider both the initial and final state radiation together for 
each scattering and impose the formation time requirement, tqcd < Af, for the inte- 
gration over the phase-space of the radiated gluons. Shuryak and Xiong [141 have 



used multiple gluon production amplitudes given by Parke and Taylor ||138|| . This 



approach has improved the gluon radiation calculation beyond the leading logarithm 
approximation. However, the space-time structure of the multiple gluon amplitudes 
is not clear and a similar LPM analysis cannot be performed. In addition, virtual 
corrections have to be included to satisfy the unitarity requirement. 

In conclusion, pQCD-based models predict the formation of a hot and dense 
parton plasma in ultrarelativistic heavy-ion collisions, in which color deconfinement 
can be achieved and gluons dominate the number of degrees of freedom. The answer 
to the question whether such a system can finally evolve into a fully thermalized 
and equilibrated parton plasma still depends crucially on the reduction of the uncer- 
tainties involved and the understanding of the underlying nonperturbative physics. 
Although theoretical progress can and must be made in these directions, the ultimate 
quantitative solution still awaits the experimental measurements at RHIC and LHC. 
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